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ABSTRACT 


This study is an attempt to optimise process conditions in plastics injection molding by 
simulation guided by genetic algorithms. The object is to finally manufacture good quality 
molded products. In the present work, design of experiments, multiple regression and the 
meta-heuristic optimisation algorithm “genetic algorithm” are used to optimise the 
simulated process. Study is motivated by the fact that the quality of an injection molded 
part aflfected by many factors that are difficult to optimise on line. These include geometric 
parameters associated with the mold design and the cooling system design as well as the 
process parameters such as the molding conditions during the filling phase. In the present 
work we have focused on molding conditions only. Quality (captured as output or 
“response” from the simulation software) is quantified as a function of flow simulation 
outputs and constitutes the objective function that must be minimised. Second-order 
regression models are developed using the simulation outputs fi'om two separate injection 
molding simulation softwares PART ADVISER (by MOLDFLOW) and 3D QUICKFILL 
(by CMOLD). The conflicting nature of the multiple objectives injection molding process 
parameters are discussed and handled by multi-objective GA. SGA and NSGA are used to 
solve the above problems. 



Chapter 1 


Introduction 


This study is an attempt to synthesise several diverse engineering analysis and optimisation 
methods to help achieve the objective of superior product quality with multiple quality criteria. 
Injection molding of plastic parts is taken as the test process. Commercial simulators are used 
to provide process performance data for stated or specified process conditions. Multi- 
objective optimisation is attempted by meta-heuristic methods and a verification of the 
integrity of the overall methodology is attempted to conclude the study. 

Chapter 2 presents an overview of multi-objective optimisation. Many real world problems 
involve multiple measures of performance, or objectives, which should be optimised 
simultaneously. In certain cases, objective functions may be optimised separately from each 
other and insight gained concerning the best that can be achieved in each performance 
dimension. However, multi-criteria solutions to the overall problem can seldom be found in this 
way. Optimal performance according to one objective, if such an optimum exists, often implies 
unacceptably low performance in one or more of the other objective dimensions, creating the 
need for a compromise to be reached (Fonesca and Fleming, 1995). A suitable solution to such 
problems involving conflicting objectives should offer “acceptable”, though possibly sub- 
optimal performance in the single-objective sense. Performance in all objective dimensions 
where “acceptable” is a problem-dependent is a subjective concept. 



The conventional solution methods that are used in solving the multiobjective problems are 
mentioned in the chapter. What makes multiple criteria decision making (MCDM) complex is 
the plurality of the criteria involved in the problem. MCDM has seemed to emerged as the 
accepted nomenclature for all models and techniques dealing with multiple objective decision 
making. This chapter concludes with the citation of the nondominated sorting genetic 
algorithm (NSGA) (Deb and Srinivas, 1995), a meta-heuristic method used in the present 
work. 

Chapter 3 describes the meta-heuristic multiobjective optimisation method based on genetic 
algorithms that is rapidly going popularity among engineers. Genetic algorithms were invented 
by John holland and were developed by Holland and his students and colleagues at the 
University of Michigan in the 1960s and the 1970s (Mitchell, 1996). This chapter provides an 
introduction to the genetic algorithms, discusses salient features of genetic algorithms and ends 
with the discussion of NSGA, the multi-objective genetic algorithm developed by Srinivas and 
Deb (1995), NSGA is used in the present work for solving a multiobjective engineering design 
problem. 

Chapter 4 overviews injection molding. Injection molding process is a complex but highly 
efficient means of producing a large variety of thermoplastic products, particularly those with 
high volume production requirements, tight tolerances and complex shapes. The final product 
quality of injection molded part must include final material properties as well as structural and 
aesthetic considerations. 

Injection molding is one of the most common and versatile operations for mass production of 
complex plastic parts with excellent dimensional tolerance. It requires minimal or no finishing 
or assembly operations. It is a cyclic process of forming plastic into a desired shape by forcing 


2 



the plastic material under pressure into a cavity. The shaping is achieved by cooling (for 
thermoplastics) or by a chemical reaction (for thermosets). Approximately 32% percent by 
weight of all plastics processed go through injection molding machines (C-MOLD Design 
guide, 1994). Historically, the major milestones in the development of injection molding are 
two: (1) the invention of the reciprocating-screw machine and (2) the application of computer 
simulation to the design and manufacture of plastic parts. 

Product quality in injection molding is a function of the material used, the mold design and the 
processing conditions. The important processing conditions during the filling phase are the 
mold temperature, melt temperature and filling time. This chapter starts also overviews the 
different aspects of mold design are discussed, and the characteristics of the plastic material 
molded. Also discussed is how process parameters are characterised. The chapter concludes 
with an outline of how NSGA may be applicable to the optimisation of injection molding. 

Chapter 5 describes the planning of injection molding simulation experiments, an effective way 
to evaluate and optimise injection molding operations before actual production begins in the 
factory. Plastics are an increasingly popular material for making parts in many industries 
including automotive, computer and consumer. The injection molding production process, a 
convenient way to make plastic parts, is nevertheless complex and extremely application- 
dependent. Previously, not all plastic parts could be analysed before production because of the 
time and expertise needed to undertake a simulation of the process before actual production 
began. This resulted in frequent and expensive part redesign and mold rework. 

But the situation is changed, it has been now recognised that computer-aided engineering 
(CAE) incorporating process simulation enhances engineers ability to handle many all aspects 
of plastic injection molding, benefiting productivity, product quality, timeliness, and cost. The 



process behaviour predicted by CAE can help novice engineers overcome the lack of previous 
experience and assist experienced engineers in pinpointing important or critical factors that 
may otherwise be overlooked. 

Keeping in view the demands of the plastics industry, giants in the injection molding software 
business, namely MOLDFLOW and CMOLD, have developed two simulation software called 
PART ADVISER and 3DQUICKFILL respectively. Both these softwares now bring the 
benefits of process simulation directly to the desk of the product designer. In the present study 
both these softwares are used for generating the data required for the optimisation of molding 
parameters. 

The techniques of multiple regression are reviewed in Chapter 6. In any system in which 
variable (input) quantities change, the interest might be in assessing the effects of the factors 
on the behaviour of some measurable quantity (the response). Such an assessment is possible 
through regression analysis. Regression analysis is a statistical technique for investigating and 
modelling the relationship between variables. Using data collected from a set of experimental 
trials, regression helps to establish empirically (by fitting some form of mathematical model) 
the type of relationship that is present between the response variable and its influencing 
factors. The ‘output’ variable is the dependent variable here and is called the response. The 
influencing factors are called predictor, regressor, or input variables. Chapter 6 also 
summarises a number of experimental plans to sample the response region. These different 
plans aim at producing good “predictor” models while keeping the experimental effort as 
economical as possible. 

Chapter 7 summarises the results of the different multi-criteria optimisation attempted in this 
study. Results obtained from executing the simple genetic algorithm (SGA) and NSGA are 



presented in this chapter. The convergence graphs from the results of SGA and the Pareto 
diagrams resulting from multicriteria optimisation by NSGA are shown. 


Chapter 8 comprises the conclusions and scope for the future work. 



Chapter 2 


MULTIOBJECTIVE OPTIMISATION 

Many real world problems involve multiple measures of performance, or objectives, which 
should be optimised simultaneously. In certain cases, objective functions may be optimised 
separately from each other and insight gained concerning the best that can be achieved in 
each performance dimension. However, suitable solutions to the overall problem can seldom 
be found in this way. Optimal performance according to one objective, if such an optimum 
exists, often implies unacceptably low performance in one or more of the other objective 
dimensions, creating the need for a compromise to be reached (Fonesca and Fleming, 1995). 
A suitable solution to such problems involving conflicting objectives should offer 
“acceptable”, though possibly sub-optimal in the single-objective sense, performance in all 
objective dimensions, where “acceptable” is a problem-dependent and ultimately subjective 
concept. In this chapter the concept of “ multiobjective optimisation” is discussed. The 
conventional solution methods that are used in solving the multiobjective problems are 
mentioned. This chapter concludes with the citation of the nondominated sorting genetic 
algorithm (NSGA) (Deb and Srinivas, 1995), which is used in the present work. 



2.1 Multiple Criteria Decision Making 


2.1.1 Definition of Multiple criteria Decision Making 

The process of decision making is the selection of an act or courses of action from among 
alternatives acts or courses of actions such that it will produce optimal results under some 
criteria of optimisation (Tabucanon, 1989). This concise definition of decision making invokes 
further elaboration to a certain extent. Before the problems can be considered well-defined, the 
set of alternatives and the set of criteria have to be known and established first; only then can the 
selection process commence. What makes multiple criteria decision making complex is the 
plurality of the criteria involved in the problem. In decision analysis of complex systems, such 
terms as “multiple criteria”, “multiple objectives”, or “multiple attributes” are used to describe 
decision situations. Often these terms are used interchangeably. Certainly there are no universal 
definitions of these terms. Multiple criteria decision making (MCDM) has seemed to emerged as 
the accepted nomenclature for all models and techniques dealing with multiple objective decision 
making or multiple attribute decision making. 

2.1.2 Examples 

Some examples of the Multiple Criteria Decision Making 

Production Planning 
Max { Total net revenue } 

Max { Minimum net revenue in any period } 

Min { Backorders } 

Min { Overtime } 

Min { Finished goods inventory } 



Capital Budgeting 
Max { Net present value } 

Min { Capital investment requirements } 

Min { Annual operating expenses } 

Max { Investment in projects related to environmental protection } 

Max { Investment in projects in a given geographical area } 

Max { Investment in projects pertaining to a given product line } 

Multiplicity of criteria occurs in almost every area of business decision making and operations. 
Examples are numerous and more straightforward than those listed so far. They appear in 
accounting, finance, operations management, marketing, manpower planning, personnel selection, 
etc., 

The general multiobjective problem requiring the optimisation of k objectives simultaneously may 
be formulated as follows. 

Max (min) Z, = ^(X), j=l,2,...,k 
Subject to 

g,(X) ^ bi, i=l,2,....,m 
X >0 

where X is a vector of decision variables and g,(X) are the inequality constraints. In multicriterion 
optimisation, because the objectives are conflicting, there does not exists a single unique solution 
which is the globally maximum or globally minimum with respect to all the conflicting objectives. 
The increase in any of these objectives will decrease the others and vice versa. 



2.2 Conflicting Criteria 


A necessary condition of MCDM is the presence of more than one criterion. The sufficient 
condition is that the criteria must be conflicting in nature. In summary, a problem can be 
considered as that of MCDM if and only if there appears at least two conflicting criteria and there 
are at least two alternatives solutions (Tabucanon, 1989). 

Criteria are said to be in conflict if the full satisfaction of one will result in impairing or precluding 
the fiill satisfaction of the other(s). The criteria are considered to be “strictly” conflicting if the 
increase in satisfaction of one results in a decrease in satisfaction of the other. The sufficient 
condition of MCDM, however, does not necessarily stipulate “strictly” conflicting criteria. 

In view of the conflicting nature of the criteria involved in MCDM, choosing the “best” 
alternatives is indeed a difficult task for the decision maker. Consequently there is a need for 
methods to systematically resolve the conflicts among criteria (or objectives) in order to reach 
acceptable compromises and come up with satisfying (or termed as “satisficing”) solutions. 

2.3 Classification of Multiobjective Problems 

There are several ways to classify the different approaches to multiobjective optimisation. 

Adulbhan and Tabucanon (1980) classified the techniques into three main approaches based on 
the way the initial multiobjective problem is transformed into a mathematically manageable 
format. These approaches are respectively, 

(a) Conversion of secondary objectives into constraints, 

(b) Development of a single combined objective function, and 

(c) Treatment of all objectives as constraints. 



Hwang, Masud, Paidy and Yoon (1982), on the other hand, proposed a different classification. 
They based their grouping of techniques according to the stage at which information from the 
decision maker is needed by the analyst. The classification is divided into four approaches , 

(a) No articulation of preference information, 

(b) “a priori” articulation of preference information, 

(c) Progressive articulation of preference information, and 

(d) “a posterior” articulation of preference information. 

Among others, who have published survey papers on multiobjective optimisation, are Johnson 
(1968), Roy (1971), Cochrane and Zeleny (1973), Lietmann and Marzollo (1975) and Hwang and 
Masud (1975). 

2.4 Solution Methods 

2.4.1 Single Objective Approach (Tabucanon, 1989) 

This is the simplest way of handling the multiobjective optimisation problems. If you consider a 
multiobjective problem with k objectives, optimise one objective (the most important one) and to 
treat the resulting k-1 objectives as “secondary” objectives into constraints. This is done by 
specifying a maximum (for minimisation) or minimum (for maximisation) level of attainment for 
each of the secondary objectives. Thus, the multiobjective problem is converted into the single 
objective optimisation problem, as follows. 

Maximise Z, = fi(X) 

Subject to 

gi(X) ^ bi, i = 1, 2, m (original constraints) 

Q(X) ^ j= 2, 3,...., k (additional constraints) 



X 


> 0 


where Zj‘ are specified minimum levels of attainment allowed of the remaining objectives. In this 
formulation it is- assumed that all the objectives are of maximising nature. Although this method is 
practical, there are certain cases where the approaches gives no defined feasible region after the 
introduction of the additional (k-1) constraints. 

2.4.2 Utility Function Method (Tabucanon, 1989) 

This utility function method converts the multiobjective optimisation problem into a single 
objective problem in the following form: 

Maximise Z = F[fi(X), f 2 (X), ..., fk(X)] 

Subject to 

gi(X) ^ 0, i = 1, 2,..., m 
X ^0 

Where F is the utility function of the multiple objectives, representing the decision maker’s 
preferences. If F is properly determined, the solutions obtained will ensure the decision maker’s 
satisfaction. However, the determination of F will be a extremely difficult sometimes. Depending 
on the problem F takes many forms, the most common form assumes that the decision maker’s 
utility function is additively separable with respect to the objectives. Thus the additive utility 
function method converts the objectives functions into one of the following form; 

Maximise Z = X 

y=i 

Fj is used in the same manner as attaching a weight to each objective function. Thus, the problem 
is transformed into the following; 



Maximise 


z = S "'.W 

M 

Where Wj indicates the relative importance of objective j, and this is determined “ a priori”. 


2.4.3 Global Criterion Method (Rao, 1994) 


In this method the optimum solution X* is found by minimising a preselected global criterion, 
F(X), such as the sum of the squares of the relative deviations of the individual objective functions 
from the feasible ideal solutions. Thus X’ is found by minimising 


F(X) 


i = l 




/,(4 


Subject to 

S(X) < 0, j=1.2,...,m 
X >0 

where is a constant and X* is the ideal solution for the /th objective function. The 
solution Xi is obtained by minimising f(X) subject to the constraints gj(X) 0,y =1, 2, ..., /w. 


2.4.4 Bounded Objective Function Method (Rao, 1994) 


In this method, the minimum and the maximum acceptable achievement levels for each objective 
function/ are specified as and l/’’, respectively, for i = 1, 2, ...., k. Then the optimum solution 
X* is found by minimising the most important objective function, say, the ^th one as follows: 
Minimise /r(X) 

Subject to 

gj(X) ^ 0, j =1, 2,..., m 

^ ,1=1,2, ...,k,i;6r 



2.4.5 Lexicographic Method (Rao, 1994) 


In this method, the objectives are ranked in order of their importance to the decision maker The 
optimum solution X is found by minimising the objective function starting with the most 
important one and proceeding according to the order of importance of the objectives The 
subscripts of the objectives indicate not only the objective function number, but also the prionties 
of the objectives Thus/i(X) and_4:(X) denote the most and least important objective functions, 
respectively The first problem is formulated as 
Minimise /i(X) 

Subject to 

g,(X) ^ 0, j =1, 2,. , m 

and its solution Xi* and fC = /i(Xi*) is obtained Then the second problem is formulated as 
Minimise / 2 (X) 

Subject to 

g,(X) ^ 0, j =1, 2, . , m 

and 

/.(X)=/i* 

Solution of this problem is obtained as X 2 * and fi ~ fiQi-i) This procedure is repeated until all 
the k objectives have been considered The /th problem is given by 
Minimise /(X) 

Subject to 

a(X) ^ 0, j =1, 2. . , m 
and 


yi(X)=/,/=l,2. ./-I 



and Its solution is found as X, and/ = /{X,*) Finally the solution obtained at the end (i e , Xt* ) 
IS taken as the desired solution X* of the original multiobjective optimisation problem 

2.4.6 Goal Programming Method (Rao, 1994) 

In the simplest version of goal programming, the analyst sets goals for each objective that he 
wishes to attain The optimum solution X* is then defined as the one which minimises the 
deviations from the set goals Thus the goal programming formulation of the multiobjective 
optimisation problem leads to 

Minimise 2 = [ (4* + djY j p ^ 1 

j=j 

Subject to 


s(X) ^ 0, 


1.2, . 

, m 

/<X) + <- 

II 

. J = 

1,2, 

d; ^ 0, 

j = 

1.2. 


d; ^ 0, 

j = 

1,2, 

,k 

d;d; = o. 

j= 

1,2, . 

,k 


where bj is the goal set by the designer for the yth objective and dj* and dj are, respectively, the 
under-achievement and over-achievement of the yth goal The value of p is based upon the utility 
function chosen by the designer. Often the goal for the jth objective, bj, is found by first solving 
the problem. 

Minimise /(X) 

Subject to 

gj(X)^0, y=l,2, .,m 


1 A 



All the conventional methods discussed above essentially convert a multiple objective optimisation 
problem into a single objective problem by some means or the other No method takes care of 
optimising all the objectives simultaneously to obtain a set of solutions (such as the “Pareto 
optimal or “efficient” solutions) which offer a way to simultaneously “satisfice” all objectives 
(Simon, 1969) Some specially constructed evolutionary algorithms (EA’s), however have been 
recognised to be possibly well-suited to multiobjective optimisation since early in their 
development Multiple individuals can search for multiple solutions in parallel, eventually taking 
advantage of any similarities available in the family of possible solutions to the problem The 
ability to handle complex problems, involving features such as discontinuities, multimodality, 
disjoint feasible spaces and noisy function evaluations, reinforces the potential effectiveness of 
EA’s in multiobjective search and optimisation, which is perhaps a problem area where 
Evolutionary Computation really distinguishes itself from its competitors (Carlos, 1995) Pareto- 
based ranking fitness assignment was first proposed by Goldberg (1989), as a means of assigning 
equal probability of reproduction to all non-dominated individuals in the population The method 
consists of assigmng rank 1 to the non-dominated individuals and removing them from contention, 
then finding a new set of non-dominated individuals, ranked 2, and so forth Srinivas and Deb 
(1995) have implemented Goldberg’s nondominated sorting in GA’s along with a niche and 
speciation method to find multiple Pareto-optimal points simultaneously, called Nondominated 
Sorting Genetic Algonthm (NSGA) This is algorithm is used in the present work for solving our 
multiobjective problem. Optimisation of Injection molding process parameter 

2.5 The Concept of Pareto Optimality 

The family of solutions of a multiobjective optimisation problem with conflicting objectives is 
composed of all those elements of the search space which are such that the components of the 



corresponding objectives vectors cannot be all simultaneously improved This is known as the 
concept of Pareto optimality 

A more formal definition of Pareto optimality is as follows consider, without loss of 
generality, the minimisation of the n components fk, k =\, , n, of a vector function /of a 

vector variable jr in a universe U is said to be Pareto-optimal if and if only if there is no jfy e U 
for which v =y(jrv) =(v;, . v„) dominates u =J{Xu) = (m/, .... mJ, i e , there is nox^e U such 
that 

V/€{1, ,n},v,^ u, n 3;e(l, ,n}|v,<M/ 

The set of all Pareto-optimal decision vectors is called the Pareto-optimal, efficient, or 
admissible set of the solutions The corresponding set of objective vectors is called the non- 
dommated set. In practice, however, it is not unusual for these terms to be used 
interchangeably to describe solutions of a multiobjective optimisation problem. 



Figure 1.1. The Efficient Front in a Bi-objective Maximisation Problem 

The notion of Pareto-optimality is only a first step towards the practical solution of a 
multiobjective problem This decision involves selecting subsequently a single compromise 



solution from the non-dormnated set according to some preference information All other 
solutions are dominated in the pareto sense and they may be ignored 



Chapter 3 


GENETIC ALGORITHMS 

Genetic algorithms were invented by John holland in the 1960s and were developed by 
Holland and his students and colleagues at the University of Michigan in the 1960s and the 
1970s (Mitchell, 1996) This chapter starts with the introduction to the genetic algorithms, 
discusses salient features of genetic algorithms and ends with the discussion of NSGA, a 
multiobjective genetic algonthm developed by Srinivas and Deb (1995), NSGA is used in the 
present work for solving a multiobjective engineering design problem 

3.1 What are Genetic Algorithms (GAs) ? 

Evolution in a changing world : Looking at the world around us, we see a staggering 
diversity of life. Millions of species, each with its own unique behaviour patterns and 
characteristics, abound Yet, all of these plants and creatures have evolved, and continue 
evolving, over millions of years They have adapted themselves to a constantly shifting and 
changing environment in order to survive. These weaker members of a species tend to die 
away, leaving the stronger and fitter to mate, create offspring and ensure the continuing 
survival of the species Their lives are dictated by the laws of natural selection and Darwinian 
evolution. And it is upon these ideas that genetic algorithms are based 

What exactly do we mean by the term GA*? Goldberg (1975) defines it as follows 

Genetic algorithms are search algorithms based on the mechanics of natural selection and 

natural genetics. 



GAs exploit Charles Darwin’s idea of the survival of the fittest and an interbreeding 
population to create a novel and innovative global search strategy A population of coded 
strings of alphabets or numbers representing solutions to a specified problem, is maintained by 
the GA The GA then iteratively creates new populations from the old by ranking the strings 
according to their fitness and interbreeding the fittest to create new stnngs, wtuch are 
(hopefully) closer to the optimum solution to the problem at hand So in each generation, the 
GA creates a set of strings from the bits and pieces of the previous strings as it occurs in the 
mating of organisms in nature occasionally adding random new data to keep the population 
from stagnating The end result is a search strategy that is tailored to probe vast, complex, 
multimodal search spaces. 

GAs are a form of randomised search, m that the way in which strings are chosen and 
combined is a stochastic process. This is a radically different approach to the problem solving 
methods used by more traditional algorithms, which often use gradient information and tend 
to be more deterministic in nature 

The idea of survival of the fittest is of great importance to genetic algorithms GAs use what is 
termed as a fitness function in order to select the fittest string that will be used to create new, 
and conceivably better, populations of strings. The fittest function takes a string and assigns a 
relative fitness value to the string The method by which it does this and the nature of the 
fitness value does not matter The only thing that the fitness function must do is to rank the 
strings in some way by producing the fitness value These values are then used to select the 
fittest strings (Holland, 1992). 

3.2 Genetic Algorithms*: A Natural Perspective 

The population can be simply viewed as a collection of interacting creatures As each 
generation of creatures comes and goes, the weaker ones tend to die away without producing 



children, while the stronger mate, combining attributes of both parents, to produce new, and 
perhaps unique children to continue the cycle Occasionally, a mutation creeps into one of the 
creatures, diversifying the population even more. Remember that in nature, a diverse 
population within a species tends to allow the species to adapt to it’s environment with more 
ease The same holds true for genetic algorithms 

Figure 3.1: The iteration loop of a Simple Genetic Algonthm 









The flowchart above shows the iterative cycle of a basic genetic algorithm. Firstly, an initial 

population of stnngs is created The process then iteratively selects individuals from the 
population that undergo some form of transformation (via the recombination step) to create 
new population The new population is then tested to see if it fulfils some stopping criteria If 
it does, then the process halts, otherwise another iteration is performed 

3.3 GAs and Robustness 

Two remarkable traits of biological systems in general are their robustness and flexibility 
Biological systems have methods for self-guidance, self-repair and reproduction Very few 
artificial systems have any of these features 

GAs demonstrate atleast some of these desirable traits from nature Atleast intuitively, we 
may think so because genetic algonthms are modelled closely on evolution occurring in the 
biological world. Genetic algorithms have proven to be robust, flexible and efficient in 
optimising vast and complex solution spaces 

3.4 Genetic Algorithm Traits 

Genetic algorithms provide robustness, efficiency and flexibility when searching a problem 
space for the optimum solution But why is this”^ For a more technical look at the power of 
GAs, a discussion on building blocks and schemata is required For the moment, we shall just 
cast a very brief look at the GA search strategy 

GAs judiciously use the idea of randomness when performing a search However, it must be 
understood that GAs are not simply random search algorithms Random search algorithms are 
inherently inefficient due to the directionless nature of their search GAs are not directionless 
They utilise knowledge from previous generations of strings in order to construct a new 



generation that will approach the optimal solution In other words, GAs use past knowledge 
to direct the search Such search algorithms are known as intelligent randomised search 
techmques 

3.5 Difference between Genetic Algorithms and Traditional Methods 

The following are the essential differences between GAs and other forms of optimisation 
(Goldberg, 1989) 

1 Genetic algorithms use a coded form of the function values (parameter set), rather than 
with the actual values themselves. For example, if we want to find the minimum of the 
function f(x) = + 5, the GA would not deal directly with x or y values, but with 

stnngs that encode these values. For this case, stnngs representing binary x values should 
be used 

2. Genetic algorithms use a set, or population, of points to conduct a search, not just a single 
point on the problem space This gives GA the power of searching noisy spaces littered 
with local optimum points. Instead of relying on a single point to search through the 
space, the GA looks at many different areas of the problem space at once, and uses of this 
information to guide it 

3. Genetic algorithms use only payoff information to guide themselves through the problem 
space Many search techniques need a variety of information to guide themselves Hill 
climbing methods require the derivative, for example The only information a GA needs is 
some measure of fitness about a point in the space (sometimes known as an objective 
function value) Once the GA knows the current measure of “goodness” about a point, it 
can use this to continue searching for the optimum. 



4 GAs are probabilistic in nature not deterministic This is a direct result of the 
randomisation techniques used by GAs 

5 GAs are inherently parallel Here lies one of the most powerful features of genetic 
algorithms GAs, by their nature, are very parallel, dealing with a large number of points 
(strings) simultaneously. Holland has estimated that a GA processing n strings at each 
generation, the GA in reality processes n^ useful substrings This becomes cleared when 
schemata are examined 

3.6 Basic Genetic algorithm operations 

The following three basic operations are found in almost every genetic algonthm 
implementation 

1 Reproduction 

2 Crossover 

3 Mutation 

3.6.1 Reproduction 

Reproduction IS usually the first operator applied on a population (Deb, 1995) Reproduction 
selects good strings in a population and forms a mating pool That is why the reproduction 
operator is sometimes known as the selection operator. The reproduction operator allows 
individual strings to be copied for possible inclusion in the next generation The chance that a 
string will be copied is based on the string’s fitness value, calculated fi-om a fitness function 
For each generation, the reproduction operator chooses strings that are placed into a mating 
pool, which is used as the basis for creating the next generation. For example, look at the table 


below 



m\fmm 

Fitness Value 


01001 

5 

19% 

1000 

12 

46% 

OHIO 

9 

35% 


From the table, it is obvious that the string 10000 is the fittest, and should be selected for 
reproduction approximately 46% of the time 01001 is the weakest, and should only be 
selected 19% of the time One always selects the fittest and discards the worst, statistically 
selecting the rest of the mating pool from the remainder of the population 


For the moment we shall look at the most commonly used reproduction method in GAs The 
Roulette wheel method simply chooses the strings in a statistical fashion based solely upon 
their relative (i e , percentage) fitness values To look abstractly at this method, consider the 
roulette wheel below, which is based on the previous example above 



Figure 3.2: Roulette wheel selection 

When selecting the three strings that will be placed in the mating pool, the roulette wheel is 
spun three times, with the results indicating the string to be placed in the pool. It is obvious 
from the above wheel that there’s good chance that string 10000 will be selected more than 
once. Multiple copies of the same string can exist in the mating pool This is even desirable. 













since the stronger strings will begin to dominate, eradicating the weaker (less fit) ones from 
the population However, there are difficulties with this, as it can lead to premature 
convergence on a local optimum 

There exist a number of reproduction operators in GA literature, but the essential idea in all of 
them is that the above-average strings are picked from the current population and their 
multiple copies are inserted in the mating pool in a probabilistic manner The following are 
five commonly used schemes 

1 Proportional (Roulette Wheel) Selection This is discussed in the reproduction section in 
detail. This method will only work with fitness values from zero (non-negative) and 
scaling may sometimes be necessary It has been shown that proportional selection 
performs poorly compared with other selection schemes in many GA problems 
2. Tournament Selection Choose t individuals at random from the population and copy the 
best individual from this group into the new population Repeat N times 

3 Truncation Selection . With truncation selection that has a threshold of T between 0 and 
1, only the fraction T best individuals can be selected They all have the same selection 
probability. 

4 Linear Ranking Selection The individual are sorted according to their fitness values and 
the rank N is assigned to the best individual, the rank 1 assigned to the worst The 
selection probability is linearly assigned to the individuals according to their rank and a 
selection equation. 

Exponential Ranking Selection This follows the same methodology of linear ranking 
selection, the only difference being that the probabilities of the rank individuals are 
exponentially weighted 



3.6.2 Crossover 


Once the mating pool is created, the next operator in the GA’s arsenal comes into play 
Crossover in biological terms refers to the blending of chromosomes from the parents to 
produce new chromosomes for the offspring The analogy carries over to crossover in GAs 

The mating pool represents the population of parent strings that survive to create child stnngs 
that survive to create child strings (Ritzel, Eheart and Ranjithan, 1994) Children are created 
in a process known as crossover Crossover begins by selecting two strings at random from 
the mating pool Next, a crossover site is randomly selected for the two strings The site can 
be between any two genes in the string The genetic material after the crossover site is then 
exchanged between the two selected strings These two new strings are considered “children” 
and form the members of a new population. The genetic operation of crossover is performed 
on each mated pair with a certain probability, referred to as crossover probability (Pc) . The 
crossover probability is typically set so that crossover probability parameter is typically set so 
that crossover is performed on most, but not all of the population The end result in the 
creation of a new population which is the same size as the old population, and which consists 
mostly of children strings whose parents no longer exists and a minonty of parents lucky 
enough to enter the new population unaltered 


100 00 
oil 10 



10010 

01100 


Figure 3.3: One point crossover 



3.6.3 Mutation 


Selection and crossover alone can obviously generate a staggering amount of differing strings. 
However, depending on the initial population chosen, there may not be enough variety of 
stnngs to ensure the GA sees the entire problem space, or the GA may find itself converging 
on strings that are not quite close to the optimum it seeds due to a bad initial population 

Some of these problems are overcome by introducing a mutation operator into the GA The 
GA has a mutation probability. Pm, which dictates the frequency at which mutation occurs 
For each string element in each string in the mating pool, the GA checks to see if it should 
perform a mutation If it should, it randomly changes the element value to new one In binary 
strings. Is are changed to Os and Os to Is For example, the GA decides to mutate bit position 
4 in the string 10000 

MUTATE 

1 0000 ► 10010 

Figure 3.4; Mutation 

The resulting string is 10010 as the fourth bit in the string is flipped 

Mutation will help in preventing the population from stagnating, mutation helps to maintain 
that diversity throughout the GA’s iterations. 

One of the major issues in GAs is the relative importance of tow genetic operators mutation 
and crossover (William, 1993) In the 1960’s, L Fogel et al. (1966) illustrated how mutation 
and selection can be used to evolve finite state automatons for a variety of tasks Rechenberg 
(1973) investigated “evolution strategies” that again concentrate on mutation as the key 
genetic operator D Fogel has continued the earlier work of L Fogel and makes an even 
stronger claim that crossover has no general advantage over mutation (Fogel & Atmar, 1990). 
On the other hand, proponents of the Holland (1975) style of genetic algorithms believe that 



crossover is the more powerful of the two operators Considerable effort has been spent in 
analysing crossover and its effects on performance (eg, Dejong 1975, Spears & Dejong, 
1991) In most of these analysis mutation is considered to be a background operator and of 
secondary importance To support these views, expenmental results have been presented, 
illustrating the power of crossover (eg, Dejong, 1975) Schaffer & Eshefman (1991) 
compared mutation and crossover, and conclude that mutation alone is always sufficient 
Upto date there has been no theoretical justification to support either groups beliefs It has 
never been theoretically shown that crossover is any sense more powerful than mutation, or 
that mutation is more powerful than crossover Similarly, no theoretical basis exists for 
supposing that both operators are necessary and perform different roles within the GA 
(William, 1993) 

3.7 Schemata : The Building Blocks of GAs 

The basic idea behind building blocks is that very fit individuals in a population pass on high 
performance notions to their children. These notions take the form of substrings It stands to 
reason that strings with a high fitness value must contain a substring that is a primary cause of 
such a high fitness Thus, even though crossover may splice the string into two, there may be 
good chance that the highly fit substring is passed on to the children. These highly fit 
substrings are known as building blocks 

Schemata are templates of strings that describe similarity between certain sets of strings In 
order to define a schema, the following alphabet is used 


Binary Alphabet 

Meaning 

0 

Binary 0 in string 

1 

Binary 1 in string 

4( 

Bits value is unspecified 




The alphabet {0, 1, * } can be used to represent any pattern of binary substrings For example 


Schema Mat 

tching strings 

* 1 1 1 1 

0 1111 

1 1 1 1 1 

*010*0 

101010 

001010 

101000 

001000 

* * 1 * 0 

00 100 


Now we look at how our basic GA operators affect the processing of schemata In 
reproduction, the effect is obvious Fitter strings get selected for the mating pool more than 
weaker strings, and thus fitter schemata have a greater chance of being involved in the 
creation of the next generation than their weaker counterparts Crossover has a huge impact 
on the GA building blocks Obviously, every time crossover occurs, there is a chance that 
useful schemata might be destroyed by the splicing process This is one of the main reasons 
why crossover should not be performed 100% of the time Looking at the following two 
schemata, we can define the length of a schema to be the distance between the first and last 
specific string position 


Schemata Length 

* ♦ 0 * ♦ * 

0 


3 

*0*1*1 

4 







From the above table, it is intuitively to see that useful shemata of long length are far more 
easily disrupted than schemata of short length 

The order of a schema to be the number of fixed positions in the schema. 


Schema order 

* 0 * * * 

1 

1 0 * 0 * * 

3 

1 1 1 ♦ * 0 

4 


The higher the order of a schema, the more specific the schema becomes obviously, “1 0 0 * 
1” is far more specific than “* * o * *” Again, with crossover, it is obvious that schemata 
with small orders have a less likely chance of being disrupted than schemata with high orders 
A high mutation rate will badly disrupt schemata, which is why a very low mutation rate is 
always advised for most GAs 

A couple of conclusions from building blocks theory are of importance to note Strings with 
higher Jitness schemata of short length will have a high likelihood of being selected to create 
the next population., and thus pass on those schemata to strings integrated the new 
population. It has been shown that schemata of this form increase in number from one 
population to the next in an exponential fashion In other words, n^ useful schemata are 
processed per generation, and the majority of these have small orders and lengths associated 
with them. These schemata are what gives a GA the power to efficiently search through a 
problem space This n^ feature is so important to GA s that it has been given a special name, 
implicit parallelism. 

The genetic algorithms implicit parallelism allows it to test and exploit large numbers of 
regions in the search space while manipulating relatively few strings. Implicit parallelism also 
helps genetic algorithms to cope with nonlinear problems, those in which the fitness of stnng 




containing two particular building blocks may be much greater (or much smaller) than the sum 
of the fitness attributable to each building block alone (Holland, 1992) 

3.8 Uses of GAs in the Real world 

Although the previous pages dealt that with GAs solely as a optimisation technique, there are 
a huge diversity of fields using GA technology for all sorts of different applications Below are 
just a small sampling of the staggering number of applications that put GAs to use 

I Criminal Suspect Recognition 
2. Music composition 

3 Construction and Training of Neural Networks 

4 Scheduling problems 

5 Games playing 

6 Prisoner’s Dilemma 

7 Earthquake Epicenter Detection 

8 Structural Optimisation 

9 Function Optimisation 

10 Database Query Optimisation 

II Aircraft Design 

12 Determination of Protein Structures 

3.9 THE NONDOMINATED SORTING GENETIC ALGORITHM : 
NSGA 

3.9.1 Introduction 

Many real-world design or decision-making problems involve simultaneous optmisation of 
multiple objectives In principle, multiobjective optimisation is very different from single- 



objective optmisation (Srinivas, 1998) In single-objective optimisation, one attempts to obtain 
the best design or decision, which is usually the global minimum or the global maximum, 
depending on' whether the optimisation problem is one of minimisation or maximisation In the 
case of multiple objectives, there may not exist one solution that is best (global minimum or 
maximum) with respect to all objectives In a typical multiobjective optimisation problem, 
there exists a set of solutions that are superior to the rest of solutions in the search space when 
all objectives are considered but are inferior to other solutions in the space in one or more 
objectives These solution are known as Pareto-optmal solutions or nondommated solutions 

Conventional optimisation techniques are difficult to extend to the true multiobjective case, 
because they were not designed with multiple solutions in mind Evolutionary algorithms, has 
been recognised to be possibly well-suited to multiobjective optmisation since early in their 
development Multiple individuals can search for multiple solutions in parallel, eventually 
taking advantage of any similarities available in the family of possible solutions to the problem 
The ability to handle complex problems, involving features such as discontinuities, 
multimodality, disjoint feasible spaces and noisy function evaluations, reinforces the potential 
effectiveness of Evolutionary algorithms in multiobjective search and optimisation, which is 
perhaps a problem area where Evolutionary Computations really distinguishes itself from its 
competitors (Fonseca and Fleming, 1995) 

3.9.2 Genetic Algorithm and Multiple Objective Optimisation 

The simplest genetic algorithm is sufficient for solving a single objective problem, with a 
single solution Multiple objective optimisation, however, entails finding a set of solutions 
along the trade-off surface between the objectives In order to accommodate a problem with 



two or more objectives, extensions to the simple genetic algorithm have been developed 
(Ritzel, Eheart and Ranjithan, 1994) 

The conceptually simplest way to use genetic algonthm to solve a multiple objective 
optimisation problem is to code the GA in such a manner that all the objectives, except for 
one, are constant (constrained to single value) The remaining objective becomes the fitness 
function for the GA Through the process of running the GA numerous times with different 
values of the constrained objectives, a trade-off surface can be developed This process can be 
tedious and time consuming In addition coding objectives as constants may be difScult or 
impossible for certain problems 

Schaffer (1984) developed a population based non-Pareto method, called Vector Evaluated 
Genetic Algorithm (VEGA) In this method appropriate fractions of the next generation, or 
sub-population, were selected from the whole of the old generation according to each of the 
objectives, separately Crossover and mutation were applied as usual after shuffling all the 
sub-populations together Although Schaffer (1984) was able to identify various solutions 
along the trade-off surface of his test functions during a VEGA a run, the algorithm always 
ultimately converged to a single solution (Ritzel, Eheart and Ranjithan, 1994) Richardson et 
al. (1989) reported that the effect of the VEGA is the same as that of using a linear 
combination of objectives in a single fitness function 

Fourman, kursaw, and Hajela and Lin, all attempted to promoted the generation of multiple 
non-dominated solutions, a goal at which they reported achieved a reasonable degree of 
success (Fonseca and Fleming, 1995) However, none makes direct use of the actual definition 
or Pareto-optimality 



Pareto-based fitness assignment was first proposed by GoIdberg(l989), as a means of 
assigning equal probability of reproduction to all non-dominated individuals in the population 
This method involves finding the set of strings in the population that are Pareto nondominated 
by the rest of the population These strings are then assigned the highest rank and eliminated 
from further contention Another set of Pareto nondominated strings are determined from the 
remaining population and are assigned the next highest rank This process continues until the 
population is suitably ranked 

Foneseca and Fleming (1995) have proposed a slightly difference scheme, whereby an 
individual’s rank corresponds to the number individuals in the current population by which it 
IS dominated Non-dominated individuals are, therefore, all assigned the same rank, while 
dominated ones are penalised according to the population density in the corresponding region 
of the trade-off surface The algorithm proceeds by sorting the population according to the 
multiobjective ranks previously determined Fitness is assigned by interpolating, e g , linearly, 
from the best to the worst individuals in the population, and then averaging it between 
individuals with the same multiobjective rank This multiobjecitive optimisation algorithm uses 
a niche-formation method to distribute the population over the Pareto-optimal region But 
instead of performing sharing on the genotype and phenotype parameter values, sharing is 
done on objective function values Hence even though this maintains diversity in the objective 
function values, this may not maintain diversity in the parameter set, a matter of importance 
for a decision maker. Moreover, this algorithm would not find multiple solutions in problems 
with different Pareto-optimal points correspond to the same objective function value 


(Jayaram, 1998) 



3.9.3 Nondominated Sorting Genetic Algorithm (NSGA) 


NSGA is based on the Goldberg’s suggestion (1989) The idea behind the nondominated 
sorting GA procedure is that a ranking selection method is used to emphasise good solutions 
and a niche method is used to maintain a stable subpopulation of good solutions NSGA 
differs from simple genetic algorithm only in the manner the selection operator works In 
NSGA Crossover and Mutation operators works same as in the simple genetic algorithm 
However, before selection is performed, the population is ranked on the basis of the 
nondominated sorting concept, to emphasis Pareto-optimality Figures 5 (Deb and Srinivas, 
1995) shows the flow chart of NSGA algorithm 

NSGA start with the generation of initial set of solutions as any GA From the initial 
population, the nondominated solutions are first identified These “nondominated” individuals 
are considered to constitute the first nondominated front in the initial population and assigned 
a dummy fitness value, equal to population size (Srinivas, 1998) The dummy fitness value is 
assigned to the initial set of population is to give equal reproductive chance to all these 
nondominated solutions in the front. To maintain diversity m the population, “sharing” (which 
is explained in the following section) is done on the classified individuals Next, the 
nondominated individuals are ignored to process the rest of the population m the same way to 
identify individuals for the second nondominated front Then dummy fitness value less than 
that of the previous front is assigned to the second front This is done to differentiate between 
the members of the first front and the members of the second front Then, sharing is done on 
the second front, and this process is continued until the whole population is evaluated, and 
classified into successive fronts. This process leads to the creation of several successive fronts 




Figure 3.5 Flowchart of NSGA 
(Deb and Snnivas, 1995) 














of “nondominatecl” individuals After the whole population has been classified crossover and 
mutation are performed on the population 

In the next generation the population has solutions according to the dummy fitness values 
assigned to the successive fronts Since individuals in the first front have the maximum fitness 
value, they always get more copies than the rest of the population, which was intended to 
search for nondormnated regions or Pareto-optimal fronts This results in quick convergence 
of the population toward nondominated regions, while sharing helps to distribute individuals 
over the nondominated region 

3.9.4 Sharing in Genetic Algorithm 

Simple GA’s have been criticised for sub-par performance on multi-modal functions They 
have been cnticised for so-called premature convergence where substantial fixation occurs at 
most bit positions before obtaining sufficiently near optimal points To overcome these 
maladies, Goldberg and Richardson (1987) has proposed a cure borrowed from nature The 
remedy causes the formation ol aiche-hke and species-like subdivision of the environment and 
population through the imposition of “sharing functions” These shaiing functions help 
mitigate unbridled head-to-head competition between widely disparate points in a search 
space This reduces the competition between distant points in a search space, which leads to 
more diverse population and more considered (and less premature) convergence 

As mentioned in the previous section, NSGA is similar to a simple GA except for the 
classification of nondominated fronts and the sharing operation A sharing function is used to 
determine the degradation of individual’s payoff due to a neighbour at some distance as 
measured in some similarity space The distance dy can be defined in two possible ways (Deb 
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and Srinivas, 1995) (1) Phenol, uc distance which is measured based on tiie difference in the 
decoded problem variables and (2) the genotypic distance which is measured based on the 
difference in the coded problem variables (in the gene space) between two individuals i and j 
The sharing of fitness in NSGA is achieved by calculating a sharing function value Sh(d,j) 
between two individuals i and j in the same front, using the following formula (Deb and 
Goldberg, 1989) 

Sh(dy) ~ 1 - (d|j / CTshare) , if dy ^ CTshare 

= 0, Otherwise 

where d,j is the phenotypic distance allowed between individuals identified by i an j, asharc is the 
maximum phenotypic distance allowed between individuals that would become member of a 
niche (There exists another type of sharing, known as genotype sharing, where the sharing 
function is evaluated using the number of different coding bits between the two individuals 
(Deb and Goldberg, 1989) Phenotype sharing has been shown to be superior to genotype 
sharing for the effective operation of NSGA (Deb and Goldberg, 1989) 

A niche (akin to biological mches formed by organisms by resource sharing) is a stable, non- 
competing sub-population of individuals surrounding the important peaks in the fitness 
landscape in the search space (Snnivas, 1998). A parameter niche count is calculated for each 
individual / by adding already calculated sharing function values (Shy) of the individual i with 
all other individuals / coexisting with it on the current nondominated front. Finally, the shared 
fitness value of each individual i is calculated by dividing its dummy fitness value by its niche 
count. 
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Chapter 4 


INJECTION MOLDING 

Injection molding process is a complex but highly efficient means of producing a large variety 
of thermoplastic products particularly those with high volume production requirements, tight 
tolerances and complex shapes The final product quality of injection molded part must include 
final material properties as well as structural and aesthetic considerations 

The Product quality in injection molding is a function of the material used, the mold design 
and the processing conditions The important processing conditions during the filling phase are 
the mold temperature, melt temperature and filling time This chapter starts with an overview 
of injection molding, in which injection molding process is discussed Different aspects of 
mold design are discussed, and the characteristics of the plastic material are mentioned Also 
discussed is how the characterisation of the process parameters is done The chapter 
concludes vrith the explanation of how NSGA can be applicable to the optimisation of 
injection molding 

4.1 INTRODUCTION 

Injection molding is one of the most common and versatile operations for mass production of 
complex plastic parts with excellent dimensional tolerance It requires minimal or no finishing 
or assembly operations. Injection molding is a cyclic process of forming plastic into a desired 
shape by forcing the plastic material under pressure into a cavity The shaping is achieved by 
cooling (for thermoplastics) or by a chemical reaction (for thermosets) Injection molding is 
one of the important processes for producing the plastic products Approximately 32% 



percent by weight of all plastics processed go through injection molding machines (CMOLD 
design guide, 1994) Histoncally, the major milestones in the development of injection 
molding include the invention of the reciprocating-screw machine and the application of 
computer simulation to the design and manufacture of plastic parts 

Since its introduction in the early 1870’s, the injection molding machine has undergone 
sigmficant modifications and improvements In particular, the invention of the reciprocating 
screw machine has revolutionized the versatility and productivity of the thermoplastic injection 
molding process Apart from obvious improvements in machine control and machine 
functions, the major development for the injection molding machine is the change from a 
plunger mechanism to a reciprocation screw The inherent advantage of using reciprocating 
machine over plunger-type machine is, it can plasticize the material more quickly and 
uniformly with its rotating motion 

4.1.1 The Injection Molding machine 

A typical injection molding machine consists of the following major components 

• Injection system 

• Hydraulic system 

• Mold system 

• Clamping system 

• Control system 

The injection system consists of a hopper, a reciprocating screw and barrel assembly, and an 
injection nozzle (Figure 4 1). This system confines and transports the plastic as it progresses 
through the feeding, compressing, degassing, melting, and injection stages Material is 
supplied to the machine in the form of small pellets or chips through the hopper throat into the 




Figure 4.1 The reciprocating-screw injection machine. 




barrel and screw assembly The barrel of injection molding machine supports the reciprocating 
plasticising screw The barrel is heated by the electric heater bands The reciprocating screw is 
used to compress, melt (by pressure) and convey the material into the mold cavity The 
reciprocating screw consists of three zones the feeding zone, the compressing (or transition 
zone), and the metering zone While outside diameter of the screw remains constant, the depth 
of the flights on the reciprocating screw decreases from the feed zone to the beginning of the 
metering zone These flights compress the material against the inside diameter of the barrel, 
which creates viscous (shear) heat. This shear heat is mainly responsible for melting the 
material The heater bands outside the barrel help maintain the material in the molten state 
The nozzle connects the barrel to the sprue bushing of the mold forms a seal between the 
barrel and the mold. The temperature of the nozzle should be set the material’s melt 
temperature or just below it 

The injection mold system consists of tie bars, stationary and moving platens, as well as 
molding plates (bases) that house the cavity, sprue and runner systems, ejector pins, and 
cooling channels (see Figure 4. 1) The injection mold is essentially a heat exchanger in which 
the molten thermoplastic solidifies to the desired shape and dimensional details defined by the 
cavity 

The hydraulic system on the injection molding machine provides the power to open and close 
the mold, build and hold the clamping tonnage, turn the reciprocating screw, drive the 
reciprocating screw, and energize ejector pins and moving mold cores A number of hydraulic 
components are required to provide this power, which include pumps, valves, hydraulic 
motors, hydraulic fittings, hydraulic tubing, and hydraulic reservoirs 



The control system provides consistency and repeatability in machine operation It monitors 
and controls the processing parameters, including the temperature, pressure, injection speed, 
screw speed and position, and hydraulic position Process control has a direct impact on the 
final part quality and the economics of the process Process control systems can range from a 
simple relay on/off control to an extremely sophisticated microprocessor-based, closed-loop 
control 

The clamping system opens and closes the mold, supports and cames the constituent parts of 
the mold, and generates sufficient force to prevent the mold from opening Clamping force 
can be generated by a mechanical (toggle) lock, hydraulic lock, or a combination of these two 
basic types 

Clamping tonnage and shot size are commonly used to quickly identify the size of the injection 
molding machine for thermoplastics Other parameters include injection rate, injection 
pressure, screw design, mold thickness, and the distance between the bars The major auxiliary 
equipment to an injection molding machine includes resin dryers, matenal-handling equipment, 
granulators, mold-temperature controllers and chillers, part-removal robots, and part-handling 
equipment. 

4.1.2 Process Cycle 

Injection molding is a cyclic process. During the injection molding process, the machine 
undertakes a sequence of operations in a cyclic fashion The basic injection molding machine 
operations, shown in Figure 4.2 (CMOLD design guide, 1994) are as follows 

• the cavity fills as the reciprocating screw moves forward, as a plunger, 

• the cavity is packed as the screw, continuously moves forward. 



Mow Closes 
Screw begins moving 
forward (or injection 




Filling comoletes 
Screw conunuii 
moving (orwara for packing 



Screw m full forward position 
Packing completes wntn gate freezes off 


G Cooling continues 

Screwtack begins to accumulate melt ahead 
of screw tip for next shot 




Cooling completes MoW 
opens for part ejection 


Figure 4.2 The basic injection molding machine operations 






• the cavity cools as the gate freezes off and the screw begins to retract to plasticize material 
for the next shot, 

• the mold opens for part ejection, and 

• the mold closes and the next cycle begins 

A process cycle is one complete operation of an injection molding machine It encompasses 
the mold closing, filling, packing, cooling, mold opemng, and ejection stages Typical process 
cycle time vanes from a few seconds to tens of seconds, depending on the part weight, part 
thickness, material properties, and the machine settings specific to a given process 

4.2 Molds 

The success of accurately producing the desired shape by injection molding process depends 
on the proper functioning of the molding machine, the processing material and ultimately the 
mold used Injection molding is usually adopted for producing thin walled components in large 
quantities at a time. So the importance of an optimal mold design is of paramount importance 
The mold design depends on the product to be made and the machine on which the mold is to 
operate. The following characteristics are to be considered in deciding the type of mold design 
for a particular injection molding process 

1 The presence or absence of external or internal undercuts on the part to be molded 

2 The type of gating and means of degating 

3 The type of ejection used 

4 The manner in which the product is released from the mold 

The various aspects that one has to taken into account in mold manufacturing are as follows 


1 Gate location 



2 Mold type selection 

3 Cavity placement 

4 Gating type 

5 Ejector system 

6 Mold material selection 

4.2.1 Gate location 

The basic principle of moldflow is that if the material flows without any hindrance to each 
extreme point in the mold then weld lines, voids and incomplete filling can be avoided In this 
context, the location of the point from which the molten metal is poured plays a very 
important role This entry point is known as the “gate”. 

Molds that can either be of the single cavity type, i e , producing one single part at a time or of 
the multi cavity type, i e., producing more than one part at a time Figure 4 3 (Crawford, 
1988) shows the possible locations of the entry point for some representative single cavity 
molds Figure 4.4 (Crawford, 1988) shows the possible locations of the entry point for some 
representative multicavity module 

The most natural and efficient entry point is the centre of the mold in case of single cavity 
molds as shown in Figure 4 3 (a) In this case, the matenal has to cover the same distance in 
all directions in the mold cavity This solution is the best even in case of multi cavity molds as 
shown in Figure 4 3 (c,d,g) For a side entry of the material, there is a danger of weld lines and 
incomplete filling shown 
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Figure 4.3 Possible gate locations for single cavity molds. 
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Figure 4.4 Possible gate locations for multi cavity molds. 




in Figure 4 3 (d,e) The possibility of incomplete filling increase in case of multi cavity molds 
If an entry point other than the center of the mold is chosen, the material cools down 
considerably by the time the most distant cavity is reached as shown in Figure 4 4 (a,b) Also, 
the pressure required for filling the mold increase substantially 

Hence, the location of the entry point or “gate” is important to have a defect free molded part 
In case of molds with simpler shapes viz rectangular or circular, the problem of gate location 
is relative easy The problem becomes complex in case of molds with shapes other than 
rectangular and circular In such cases, the necessity of an optimal location of the gate arises 
The location must be based on the geometrical features of the mold 

4.2.2 Mold type selection 

The type of mold to be used depends on the geometrical attnbutes viz shape and size of the 
final product The type of mold may either be of the two plate mold type or the three plate 
mold type 

The two plate type is the simplest type of mold, comprising essentially of two die plates which 
carry the female cavity and the male punch respectively. Such types of molds are used for box 
or bowl type moldings. The three plate mold type has a third plate, generally known as the 
floating or centre plate in addition to the two plates as in two plate mold type Three plate 
mold types are frequently used for the production of small or medium sized components 
Figure 4 5 and Figure 4 6 show the two plate mold type and three plate mold type 
respectively 

Depending on the geometry of the part, the location of the gate and the machine on which it is 
to be used, a suitable selection of the type of mold is done 




Figure 4.5 A two plate mold 
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4.2.3 Cavity placement and Gating type selection 


The basic reason for the employment of multiple cavities is to obtain a greater output in a 

given time The placement of these cavities has a considerable effect on the final product 

characteristics The prominent types of cavity placements are star, asymmetrical and 

symmetrical Depending on the location of the gate and the machine parameters of clamping 

force and platen size, a proper arrangement of cavities is done Figure 4 7 shows the layouts 

possible with the associated runner system for injection molds 
Molded pari 
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Symmelrical circular 
cavity layout 



Molded part 



Asymmetrical 
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Figure 4.7 Cavity Layouts in injection molds 
The type of gating depends on the mold geometiy and the use to which the final product might 
be put to The selection of the gating system is to be done from the following types of gates 
direct, pinpoint, film, rectangular 



4.2.4 Ejector System Selection and Mold Material Selection 


Injection molded compounds exhibit a certain amount of shrinkage, i.e , in a cooled condition 

the volume is somewhat smaller than m a heated condition If parts with accurate dimensions 

are needed, allowances must be made for this shrinkage when establishing the dimensions for 

the cavities Shrinkage also causes the molded part to sit tight on the cores The parts are 

ejected after they have cooled down As a result, special ejection methods like pins, stripper 

nngs, air and slides are to be appropriately matched with the part geometiy and machine on 

which the ejector system works CENTRAL library 
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Various materials have been reporteJ^o^ used ^r^ol^ manufacturing and their selection 
mainly depends on the process of manufacturing the mold, i e , whether a conventional method 
like milling is used or a conventional method like EDM is used A suitable selection can be 
made depending on whether the material in question is to be used to make the base, inserts, 
plates or cores The usual materials are heat treated steels, plain tool steels, nitriding steels 

4.3 Plastics 

Plastics are one group of polymers which are built from relatively simple units called 
monomers (or mers) through a chemical polymerisation process Processing polymers into end 
products mainly involves a physical change (for thermoplastics) or a chemical reaction (for 
thermosets) ThefollowingASTM definition of plastics is widely accepted (Watson, 1994). 

A material that contains as an essential ingredient one or more 
polymeric substances of high molecule weight. It is solid in its 
finished state and, at some stage m its manufacture or 
processing, can be shaped by flow into finished articles. Synthetic 



fibres, textiles, natural and synthetic rubbers, adhesives and 
paints which may m some cases comply with this definition, are 
not considered to be plastics 

The basic structure of a polymer molecule may be visualized as a long chain of repeating units, 
with additional chemical groups forming pendant branches along the primary “backbone” of 
the molecule The structure arrangement, size, and chemical constitution of the polymer 
molecule have a direct influence on its physical and chemical properties In addition, the 
macromolecular nature of the plastics implies that their material properties may also be 
dependent on the mechanical and thermal histoiy that the materials expenence during 
processing For example, the viscosity (which indicates the material’s resistance to flow) of a 
polymer melt increases with increasing molecular weigh, but decreases as temperature 
increases Further, the aligned molecular orientation that results from strong shear exerted on 
the material also reduces the viscosity of the polymer melt. The physical and mechanical 
properties, as well as the cost of polymers, can be modified by blending a number of polymers 
or by compounding polymers with other materials or reinforcing agents 

4.3.1 Classification of Plastics 

Based on the type of chemical reaction (polymerisation) that links the molecules together, 
plastics are classified as either thermoplastics or thermosets 

4.3. 1.1 Thermoplastics 

Thermoplastics typically have high molecular weights resulting from a high degree of 
polymerisation The long molecular chain, either linear or branched, has side chains or groups 
that are not attached to other polymer molecules. As a result, thermoplastics can be repeatedly 



softened (or hardened) by an increase (or decrease) in temperature This type of phase change 
without a chemical reaction permits the recycling of thermoplastic scraps such as the trimmed- 
ofF runners and sprues from injection molding Although thermoplastics are recyclable, it is 
very likely that a small degree of chemical change (e g , oxidation, thermal degradation) takes 
place during processing, and therefore the properties of recycled polymers may not be 
equivalent to those of the virgin polymer 

Thermoplastics account for more than 70% of all polymers produced Among thermoplastics, 
the commodity resins, are high-density polyethylene (HDPE), low-density polyethylene 
(LDPE), polypropylene (PP), polystyrene (PS), and polyvinyl chlonde (PVC) These materials 
constitute 90% of all thermoplastics 

4.3.2 Thermosets 

Cross-linking is a chemical process in which chemical bonds form among molecules of 
thermosetting materials, resulting in an interconnected network The cross-linking process is 
the principle difference between thermoplastics and thermosets Pnor to molding, the chain- 
hke structure of thermosets are similar to thermoplastics During processing, thermosets 
polymerize (react or cure) with the activation of heat and/or a chemical means into a cross- 
linked microstructure Once the reaction is completed, the polymer chains are bonded (cross- 
linked) together to form a three-dimensional network These cross bonds among molecules 
prohibit the slippage of individual molecular chains. Consequently, a thermoset becomes an 
infusible and insoluble solid and cannot be re-softened and reprocessed through the application 
of heat, without degrading some linkages. Typical examples of thermosets are, pheno- 
formaldehyde, melamine-formaldehyde, urea-formaldehyde and epoxies 



Table 4.1 Comparison of structures and properties of thermoplastics and thermosets 


Material 

Thermoplastics 

Thermosets 

Microstructure 

Linear or branch molecules 

No chemical bonds among the 

molecules 

Cross-linking network with 

chemical bonds among 

molecules after the chemical 

reaction 

Reaction of Heat 

Can be re-softened (physical 
phase change) 

Cannot be re-softened after 

cross-linking without 

degradation 

General Properties 

Higher impact strength. 

Easier processing 

Better adaptability to complex 

designs 

Greater mechamcal strength 

Greater dimensional stability 

Better heat and moisture 

resistance. 


4.4 Processing Parameters 

The mechanical and physical properties of a molded part do not depend solely on the chemical 
constitution of its material The processing conditions exercise a considerable influence on the 
strength, distortion, dimensional stability, surface finish and tendency to stress cracking These 
factors that represent the part quality are not apparent externally, but are reflected by the 
structure of the part. 

The variables that affect the Injection molding process can be categorised into geometrical and 
operational considerations Geometrical considerations such as gate location affects the part 
quality by affecting the temperature, pressure and shear stress by alternating the balance and 


direction of the material flow 










Operational considerations such as melt temperature, mold temperature and fill time affect part 
quality by changing the conditions under which the mold is filled The operational 
considerations are said to make up the Injection molding conditions Also the selection of an 
appropnate cooling system is of immediate importance since an improperly selected cooling 
system can result in hidden stresses in small parts, distortion in larger parts with relatively thin 
walls and even crack formation Hence, an approach to determine the suitable cooling system 
and the optimal combination of temperature, pressure and time is needed 

4,4.1 Characterisation of the Process 

A number of process parameters may affect the formation of a good molded part However 
the major ones are 

• Pressure 

• Temperature 

• Time 

In the context of characterisation of the process, the parameter pressure is usually referred to 
as the cavity pressure The distribution of the cavity pressure at different stages of the molding 
process as it proceeds with time immensely influences the quality of the molded part Figure 
4 8 (Johnnaseri, 1983) shows the cavity pressure profile over time 

The information available from the cavity pressure curve can be divided into three stages. 

1 Filling of the cavity ( Injection stage) 

2 Compression of the melt ( Compression stage) 

3 Holding the solidifying material under pressure (Holding stage) 
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Figure 4.8 Cavity pressure profile over time 

These three stages can be related to quality criteria effects. The injection pressure has the 
function of overcoming the flow resistance from the nozzle to the cavity The compression 
pressure plays a significant role in avoiding defects like shrinkage and voids The holding 
pressure is responsible for maintaining an uniform pressure However, as the mold is 
solidifying, the pressure decreases The major effects of an improperly maintained holding 
pressure are voids, sink marks and dimensional instability 

The most important processing conditions during the filling phase are the mold temperature, 
the melt temperature, and the filling time ( Pandelidis and Zou, Part-I, 1990) The selection of 
correct molding conditions in the injection molding process is a non-trivial task An increase in 
melt temperature causes a decrease in melt viscosity, which results in reduce pressure 
requirements and reduced stresses. On the other hand, high melt temperature may increase the 
possibility of material burning and will also increase cooling time Increasing the mold 






temperature reduces heat losses, and the maximum temperature difference at the end of the 
filling phase may be reduced However, a high mold temperature increase cooling time. Short 
fill time require higher pressure because a higher flow rate requirement, resulting in higher 
shear rate and shear stress Then, if the fill times are too long, pressure will increase because 
the plastic temperature decreases and the viscosity increases It is clear from above discussion 
that some optimisation to balance the conflicting processing parameters is required to produce 
a good quality product (Pandelidis and Zou, Part-II, 1990) 

The time for one single run usually is called the cycle time The cycle time consists of filling 
time, dead time and cooling time As cycle time is dependent on the filling time and rate of 
cooling, faster rate of filling and cooling will results in increase of production rate For this, 
the cooling system must provide the most umform possible temperature distribution in the 
mold so that cooling is uniform and homogeneity in microstructure results 
4.5 Quality Measures and Processing parameters optimisation 

It is required to develop quantitative measures of the part quality since the ultimate goal is to 
improve the final product characteristics The term quality may be referred to many product 
properties such as mechamcal, electrical, optical or geometrical There are two types of part 
quality measures direct and indirect (Pandelidis and Zou, Part-I, 1990) 

A model that predicts the warpage from flow simulation results would be characterised as a 
direct quality measure. In contrast, an indirect measure of quality is a quantity that is co- 
related but does not produce a direct estimate of that quality However, direct quality 
measures are usually more expensive to evaluate than indirect quality measures Prediction of 
warpage is too time consuming. So, an alternative measure is adopted 



Many factors affect the quality of the final part One of the most damaging causes is due to 
warpage Warpage in the part occurs due to 

1 Uneven temperature distribution and uneven cooling in the part 

2 Inhomogeneity of part density 

3 Anisotropy in thermal shrinkage 

4 High unbalanced residual stress distribution 

The major causes of material degradation are high temperature (burning) and high shear rate 
The indirect quality measures are those related to warpage and material degradation Instead 
of using finite element stress and thermal analysis to predict warpage and material degradation 
directly, we quantify the major causes of warpage and material degradation as described in the 
subsequent sections For optimisation, we restrained ourselves with quality measures which 
are controllable at the filling phase The packing and cooling phase must also be considered 
when full design is implemented 

4.5.1 Temperature Difference 

Temperature difference is a measure of the temperature distnbution uniformity in the mold 
once it has been filled by the material The temperature used here is the average temperature 
The amount of shrinkage that a plastic undergoes during cooling is a function of its 
temperature (Pandelidis and Zou, Part-I, 1990) All other factors being equal, a nonuniform 
temperature distribution will cause differential shrinkage and, therefore, warpage in the 
cooling process Ideally, the temperature at each point should be same at the end of filling 
This ensures the non uniform shrinkage due to difference in temperature will be eliminated 
Unfortunately, this rarely occurs. Heating occurs in the flow because of fnctional and shearing 
effects At the same time, cooling takes place at the plastic mold interface. Thus, in general, at 



the end of fill stage there will be a nonisothermal temperature distribution throughout the part 
In a finite difference simulation, there will be discrete temperature differential The subsequent 
differential shrinkage that occurs \sall induce residual stresses causing warpage and other 
detrimental short and long-term effects 

The parameter Tj has been defined to quantify the severity of the temperature difference Td is 
simply the difference between the maximum and minimum temperature at the end of fill over 
entire part which are available in simulation results This provides a reliable measure of 
uniformity of temperature and can be represented mathematically as follows 

Td = Tn^.-Tnvm (4 1) 

where, 

Tmax = maximum temperature in the mold 
Tmin = minimum temperature in the mold, 

4.5.2 Pressure 

Another cause of warpage is the existence of non-uniform pressure in the mold after it has 
been filled. During the cooling phase the areas with relatively high pressure will shrink less 
than areas of lower pressure 

There will be pressure losses dunng the filling stage These losses are the result of the 
continual shearing and fnction that occurs as the polymer melt is forced into the mold Ideal 
conditions will exists if in each part of the mold has the same pressure gradient Higher 
pressure will result in higher shear rate and shear stresses which will cause warpage. 
mechanical sticking and flash (John Bown, 1 982) 

The parameter P^ax has defined to quantify the severity of pressure distnbution 



4.5.3 Optimisation 


Process Variable 

Effect 

Quality Problems 

Melt temperature 

Increasing the melt temperature causes a 
decrease in melt viscosity, which results in 
reduce pressure requirements and reduced 
stresses On the other hand, high melt 
temperature may increase the possibility of 
material burning and will also increase cooling 
time 

Matenal degradation, 
and Short shot 

Mold temperature 

Increasing the mold temperature reduces heat 
losses, and the maximum temperature 
difference at the end of filling phase may be 
reduced. High mold temperature is required for 
improved part appearance is a concern. A 
higher mold temperature produces higher gloss 
and more crystallization However, a high mold 
temperature increase cooling time 

But if rapid curing takes 
place the moldings are 
liable to exhibit porosity in 
the thicker sections and 
some burning may takes 
place Flow marks, 

warping and excessive 
shrinkage may occur with 
too rapid curing 

Fill time 

Short fill time require higher pressure because a 
higher flow rate requirement, resulting in higher 
shear rate and shear stress Then, if the fill 
times too long, pressure will increase because 
the plastic temperature decreases and the 
viscosity increases 

Warpage, Excessive 

shrinkage and flashing 


From the discussion above we can say that, the quality of the injection molding part is affected 
by the various factors of mold design, filling phase, cooling phase and packing phase. In this 
study we are considering the parameters related to the filling phase i e , processing conditions 
that effect the quality of the injection molding part From the discussion above it is clear that 
the process parameter optimisation of the injection molding part is a multiobjective 
optimisation problem. As the objectives are conflicting (discussed in section 4 4 1), we can 
apply NSGA for solving the problem 








Chapter 5 


MOLDING SIMULATION EXPERIMENTS 

Plastics are an increasingly popular material for making parts in many industries including 
automotive, computer and consumer The injection molding production process, a convenient 
way to make plastic parts, is nevertheless complex and extiemely application-dependent. 
Previously, not all plastic parts could be analysed before production because of the time and 
expertise needed to undertake a simulation of the process before actual production began This 
resulted in frequent and expensive part redesign and mold rework 

But the situation is changed, it has been now recognised that computer-aided engineering 
(CAE) including pi ocess simulation enhances engineers ability to handle all aspects of plastic 
injection molding process, benefiting productivity, product quality, timeliness, and cost. The 
process behaviour predicted by CAE can help novice engineers overcome the lack of previous 
experience and assist experienced engineers in pinpointing factors that may otherwise be 
overlooked. 

Keeping in view the demands of the plastics industry people to giants in the injection molding 
software business, namely MOLDFLOW and CMOLD, have developed two simulation 
software called PART ADVISER and 3DQUICKF1LL respectively. Both these software’s 
now bring the benefits of process simulation directly to the desk of the product designer. In 
the present study both these software’s are used for generating the data required for the 
optimisation of process parameters of injection molding Discussed in Chapter 4. 



Both the software’s accept solid geometry models from CAD software, built from any of the 
popular CAD systems such as Pro/Engineer, IDEAS, or Unigraphics etc., CMOLD’s 
3DQU1CKFILL lakes the injection molded pait design oftlic stl foiniat whereas Partadviser 
accepts the part in mpa format 

5.1 C-MOLD 3D QUICKFILL 

5.1.1 Introduction 

C-MOLD 3D Quickfill (Quickfill) is a desktop plastics CAE software tool that helps part 
dcsigneis, mold designeis, and moldeis incoiporate engineering decisions and 
manufacturability considerations into the earliest stages of product development Vendors of 
CMOLD claim that the mfoimation provided by Quickfill helps designers test feasibility, 
manufacturability, cost, and performance concerns against real-life constraints As mentioned, 
C-MOLD 3D Quickfill accepts solid geometry models from CAD software’s like Pro- 
engineer, I-DEAS, Unigiaphics etc , in stl (steieolithography) format. 

5.1.2 Procedure to start process simulation on the C-MOLD 3D Quicltfill 
Software 

The following are the steps that one has to follow to run the simulation on Quickfill 
1. Start the Quickfill : Quickfill runs both in the Windows and Unix environments 

2 Open your solid model : Quickfill accepts solid geometry from any CAD product, saved 
in the STL ( stereolithography) format One can import an STL file or can open a surface 
mesh in CMOLD FEM format 

3 Control the view : Quickfill provides the controls so you may position the model for the 
best view, using the vai ious control buttons. 



5.1.3 Features of CMOLD 3D Quicldill 


1 For the process conditions used for the simulation, Quickfill gives whether the part either 
short or it doesn’t and Quickfill will show you exactly where the short shot occurred. If 
your part short shots, Quickfill tells you the remedy options If your design requires 
excessive injection picssuie or melt temperature, it tells what can be done. It even tells you 
when, and by how much, the part’s thickness can be reduced and still be manufactureable. 

2 Quickfill supplies the simulation results vital to knowing why some problem exists. 
Simulation results show how the four variables of injection molding (material, part, mold 
and machine) influence part quality 

3 For parts with multiple gates, mold designers would like to know how good each gate is 
performing Quickfill tells you the % cavity filled through each gate This allows you to 
strategically remove the under-performing gates and assist over-performing ones. 

4 Using the simulation output of molecular and fiber orientation you can significantly 
improve part strength in ciitical areas (especially in the vicinity of weldlines). Quickfill 
helps you design for optimal performance. Part warpage can be minimised by selecting 
gate locations which balance the filling of the mold cavity. 

5. In the Design advice given by Quickfill after each simulation, includes the injection 
pressure and machine clamp force tonnage required to mold the part, as well as 
suggestions for correcting or improving the design. For example, Quickfill tell you whether 
part thickness can be i educed, and by loughly what percentage it can be reduced, while 
still producing a part that will fill Quickfill will also tell you if the variation in part 
thickness is large, oi if it is too huge 



6 Quickfill shows the design specifications used for the simulation The specifications and 
calculated icsults include 

• the resin selection and the Melt mass-flow late (MFR) 

MFR IS a world-wide industry standard that can be thought of as a measure of how easily the 
melt can flow MFR, is an ISO standaid, is identical to the melt-flow index (MFI), an ASTM 
standard 

• the machine clamp force and shot size required for the part 

Shot size is a mcasiiie of an injection molding machine’s capacity. The measurement is the 
product of the area of the screw (ram) and the maximum travel that it can achieve. 

• machine setup information, including the filling and cooling time, the required injection 
pressure and the melt and mold temperatures 

• Part information including the part weight and projected area of a single cavity, the clamp 
force required for single cavity, the number of cavities, the parting plane, and the averaged 
part thickness and standard deviation of the part thickness, If the standard deviation is 
greater than 5%, the Design Advice report will warn you that the thickness variation is large; 
if It IS greater than 10%, you will be warned that the thickness variation is too large. 

5.1.4 SIMULATION RUNS 

In this section simulation results and some of the plots obtained from the Quickfill are 
presented with the help of one trial run. For showing the results obtained from Quickfill, a 
“tape” model is imported from the inbuilt database. LDPE plastic resin is chosen as the for this 
trial run Default process conditions are used for the simulation run although if we want, we 
can change the process condition values such as mold temperature, melt temperature, number 



of cavities and machine builder values of maximum injection pressure and maximum injection 


rate 


Figure 5 1 shows the melt front advancement plot, this plot shows the filling pattern and 


whether the fill is balanced Fill time is noted from this plot 
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Figure 5 3 shows the temperature plot Temperature plot shows by degree of greyness the 
spatial temperature variation across the part at the end of filling Greater the temperature 
variation on a part, the more potential there is for the part to warp From this plot temperature 
difference is noted. 



Figure 5,3 


Temperature plot 







Figure 5.5 shows the Dc.„gn ajv.cu given 6y the QuickOII for the above trial run. The advice 
includes the injection pressure and machine damp force tonnage require to mold the part, as 
well as suggestions for correcting or improving the design 


Date . NOV03-98 
Time ; 07:37:34 
Name . tape 


C-MOLD 3D QuickFill 98.6 Design Advice. 


The part(s) can be successfully filled with injection pressure 
of 32 MPa and clamp force tonnage of 12 ton(m). 

The clamp force tonnage does not include the sprue/runner/gate systems, 
inclusion of the sprue/runner/gale system will increase the required 
clamp force. m/ 

No resin degradation detected in the simulation. 

The required injection pressure is in a reasonable range, or 18 % 

of the max injection pressure which can be delivered by a typical machine. 

The part design has very uniform thickness. 

The part thickness can be reduced by roughly 40 % and the part will still 
fill. Before reducing the thickness, however, you should consider the 
minimum thickness required for part strength. If you reduce the thickness, 
rerun 3D QuickFill for updated results 


Design Advice 


l‘"igure 5.5 



Figure 5.6 shows the Design specifications given by the Quickfill for the trial run. 


Design Specifications 


Date NOV03-98 
Time 07 37 34 
Name tape 

==== Resin specification ==== 

Polymer PE-LD (base) 

==== Machine specification ==== 

Required clamp force = 12.33 ton(m) 

Required shot size (GPS) = 11.48g 

The clamp force and shot size estimates do not include Ihesprue/runner/gate system 


Inclusion of the sprue/runner/gate system 
shot 
size 

Max injection pressure 
Max injection rate 

Machine performance = 

==== Machine setup ==== 

Fill time = 

Cooling time = 1.50 

Injection pressure 
Melt temperature 
Mold temperature 

==== Part information ==== 

Single cavity part weight 
Average part thickness 
Standard deviation of part thickness 
Single cavity projected area 
Required clamp force for single cavity 
Number of cavities 
Parting plane is in XY plane 
1000% filled by gate# 1 


increase the required clamp force and 


= 180.00 MPa 

= 26 61 cm'^S/s 

Good 


0 44 s 

= 31.64 MPa 

= 220 00 C 

= 40.00 C 


9 89g 
1 19 mm 
0.01 mm 
47.79 cm''2 
12 33 ton(m) 
1 


Figure 5.6 Design Specifications 



5.2 MOLDFLOW - PART ADVISER(3.0) 


5.2.1 Introduction 

Part Adviser is a simulation software of the Moldflow Corporation, one of the leaders in the 
process wide plastic simulation This software, the vendor claims, enables the quick evaluation 
of each and every injection molded part design for “manfacturability” 

The Part Adviser (PA) reads the model of your part created by a standard CAD package PA 
analyses a CAD solid model directly and eliminates the need to translate geometry using an 
intermediate file or perform complex FEA meshing Some inputs are required, including the 
material to be used, and the injection location 

5.2.2 Procedure to be followed to start a molding process simulation on PA 

The following steps are to be followed to run the simulation on a PA . 

1 Open the PA from the desktop : Open the PA after double clicking the icon from the 
desktop (PA is compatible to Windows95 and Windows NT 4.0) 

2 Open a model PA accepts the solid models of the part made from the CAD packages 
such as I-DEAS, Pio-Engineer and Unigraphics etc.. Open the model from stored 
database 

3 Select a polymer All data related to polymer materials are stored in two locations, the 
standard database and local database. If an exact match of supplier and tradename cannot 
be found in the database, alternative materials may be selected using the polymer selection 
option and match criteiia Four criteria are used to search the database. They are listed in 
order of priority as below 


polymer family 



9 


melt ilow index 


• type of filler 

• percentage of filler 

If you provide the search data, the matches found and a window is dynamically updated, 
showing mateiials that match the criteria provided 

4. Manipulate the model . There are number of tools available to rotate and manipulate the 
model in PA These can assists in detei mining a best location for polymer injection to view 
results Various tools available are “Pan the model”, Rotate the model”, Zoom in model”, 
“Banding zoom”, “Select a view” and “fit the model view on the window” 

5 Selecting an polymer injection location The PA allows to work with single or multiple 
polymer injection locations 

6 Changing analysis results ; The PA allows you to change some of the model processing 
conditions In the “material properties” area, we can change the mold temperature and the 
melt temperature fiom the defaults to custom values In the maximum pressure injection 
limit area, we can change the maximum pressure injection limit from the default values to a 
custom value, and we can also change the default value. In the “machine injection time 
area” we can entci a paiticular time value. These changes are not compulsory steps for 
each analysis 

7 Save the model Saving will ensure that all the conditions you have set so far will be 
saved They includes polymer selected, the rotation of the model, the injection location and 
the current maximum injection pressure limit 

8 Run the PA In this step we run the PA to fill the model The filling process will display 
dynamically on the 3d model and show the plastic flow Filling the model will take 
appioximately 5-10 minutes Duimg filling, the status bar at the base of the screen will 



display the cm i cut i^iogicss bar (o iiidicalc progicss. When filling has been completed the 
status bar will display ready 

9 Results : Various results that can be viewed after the analysis 

• Confidence of fill 

• Fill Time 

« Pressure Drop 

• Flow front Temperature 

• Weld lines 

• All tiaps 

• Injection Piessuie 

5.2.3 Features of Part Adviser 3.0 

1 After the analysis has finished, you can look at the confidence of fill result, which is 
derived fiom the filling time, injection pressure, pressure drop and temperature. The 
confidence result is displayed in thiee colours. A green area has a “high” confidence rating, 
yellow indicates a “medium” confidence rating, while red shows a ”low” confidence rating 
Theie aie nuniboi of confidence of fill advice help topics that give suggestions on how to 
Improve a medium oi low result If the confidence of fill result does indicate problems, we 
can look at the filling time, injection pressure, pressure drop and temperature results to 
determine the reason for the low or medium confidence of fill If the confidence of fill 
results shows any molding problem, we can identify the problem after looking at the 
pressure and the temperature results We can find out the other problems after looking at 
the Weld lines and An trap results 



2 PA provides advice on specific aspects of the filling process, provide explanations of 
potential problems allow you to make certain changes that may lead to solution for filling 
difficulties 

3 For clarity of the users PA has clipping plane to view “hidden” sections of the model as 
well as an X-Y-Z rotation axes display 

4 PA will provide with the enhanced web based reports after the analysis is completed 

5.2.4 Simulation Results 

In this section various lesults and some of the plots obtained by the PA is presented with the 
help of a trial run As mentioned earlier PA generates a report after each simulation is 
completed, which contains the results summary and the various plots like temperature plot, 
pressure diop plot etc, to interpret the results. For this section we have taken the example 
“Diskeady” as the part and the various settings are made as explained in the section 5.2.2 and 
the simulation is done 

Figure 5 7 is the injection pressure plot given by the PA after the process simulation. In this 
plot colour (which can not be viewed in this black and white picture) at each place on the 
model represents the present at the place on the model, at the moment of the part is filled 
completely This is a “snapshot” result, that is it shows the pressure through the whole part at 
the end of fill In this plot, range of colours are used to indicate the region of lower pressure 
(coloured blue) through to the region of highest pressure (coloured red) From this plot 
maximum injection pressure (PmaO is noted 


Figure 5 8 shows the changes in the temperature of the flow front during filling. This result 
uses a range of colours to indicate the region of lowest temperature (coloured blue) through to 
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the region of highest temperature (coloured red). From this plot temperature difference is 
noted, by taking the difference between the maximum and minimum temperatures (indicated in 
the scale situated on the right side of the plot). 


Figure 5.8 Flow front temperature plot 
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In jection pressure plot 





Figure 5 9 IS the fill time plot This plot shows the flow path of the plastic through the part. 
Each colour shown (which can not be viewed in this black and white picture) represents the 
parts of the mold which wcie being filled at the same time. A short sliot (a part of the model 
which did not fill) will be displayed as tianslucent. In this plot range of colours used to indicate 
the first region to fill (coloured red) through to the last region to fill (coloured blue) The 
value of actual fill time is noted down from this plot 



Figure 5.9 Fill time plot 



CHAPTER 6 


MULTIPLE REGRESSION 

6.1 Introduction to Regression analysis 

In any system in which variable quantities change, the interest might be in assessing the effects 
of the factors on the behaviour of some measurable quantity (the response). Such an 
assessment is possible through regression analysis Regression analysis is a statistical technique 
for investigating and modeling the relationship between variables. Using data collected from a 
set of experimental tiials, regression helps to establish empirically (by fitting some form of 
mathematical model) the type of relationship that is present between the response variable and 
its influencing factors Tlie response variable is the dependent variable and is called the 
response, and the levels of the influencing factors are called predictor, regressor, or input 
variables Regression analysis is one of the most widely used tools for investigating cause-and- 
effect relationships having applications in the physical, biological, and social sciences, as well 
as m engineering and in many other fields (Montogomery and Peck, 1992) 

By predictor oi independent vaiiablcs we shall usually mean variables that can cither be set to 
be a desired value or else take value that can be observed but not controlled. Independent 
variables in a regression analysis can be qualitative as well as quantitative (Khuri and Cornell, 
1987) The specific independent variables whose levels are to be studied in detail are those 
that are quantitative in nature, and their levels (or settings) are assumed to be controlled 
(without error) by the experimenter As a result of changes that are deliberately made, or 
simply take place in the predictor variables, an effect is transmitted to other variables, the 
} espouse variable 



The distinction l)ctwccn picdiclui and response vaiiables is not always completely clear cut 
and depends sometimes on oiii objectives. What may be consideicd a i espouse variable at a 
midstage of a process may also be regarded as predictor variable (Draper and Smith, 1981). 

6.2 Multiple Regression analysis 

A regression model that involves more than one regressor (independent) variable is called a 
multiple regression model. When there are more than one independent variables influencing 
the response variables, the method adopted to build the model is the multiple regression. 
Before one starts to build an empiiical model, it is very necessary to understand the range (of 
variables) in which model hijs to be developed, because developed model may be linear in 
some part of the range and nonlinear in some other part of the range. However if range is very 
large there is less chance for the model to be linear In real life situations it is very rare that one 
encounters a linear model, but at the same time, a second order model is usually good enough 
to capture a great deal of information, To build a second order model, at least a three-level 
experimental design or plan is needed For a five-factor experiment, 3^ experiments are needed 
for Si full-faclonal three level design (Jugal, 1997). Thus full factorial design will certainly 
leads to a veiy accurate model (Jugal, 1997), but at the same time it will often be costly to do 
a such a large number of experiments. So, it is always the desire of the experimenter to choose 
a design which was fewer number of experimental runs, while at the same time, a high quality 
of information is generated from such an experimental design. Fractional factorial designs are 
therefore used to carry out the model estimation experiments This does reduce the number of 
experimental runs, while giving good deal of good infoimation about the process (Box and 
Behenkin, 1960) 



Various designs used Ibi developing second order models are described in section 4.4. 

6.3 Response Surface Methodology (RSM) 

RSM is a collection of mathematical and statistical techniques that are useful for modelling and 
analysis in applications where a response of interest is influenced by several variables and the 
objective is to optimise the response The general RSM approach was developed in the early 
1950’s, and from last twenty years RSM has found extensive applications in a wide variety of 
industrial settings, ranging from chemical processes, semiconductor and electronic 
manufacturing, machining, metal cutting, and joining processes, among many others (3 1). 

RSM IS a set of techniques that encompasses (Khun and Cornell, 1987) 

1. Setting up a series of experiments (designing a set of experiments) that will yield adequate 
and reliable measurements of the response of interest 

2. Determining a mathematical model that best fits the data collected from the design chosen 
in (1), by conducting appropriate tests of hypotheses concerning the model’s parameters 
and 

3 Determining the optimal settings of the experimental factors that produce the maximum 
(or minimum) value of the response 

If discovering the best value, or values, of the response is beyond the available resources 
of the experiment, then response surface methods are aimed at obtaining at least a better 
understanding of the overall system (Khuri and Cornell, 1987), 

Box and Draper (1975) have listed 14 properties of a response surface design to be used when 
fitting a polynomial model to data collected at the design points. Some of the features of a 
desirable design are as follows (Montogomery, 1997), Such a design 



1 Provides a reasonable distiibution of sampling or data points ( and hence information) 
through out the region of interest 

2 Allows the model adequacy, including lack of fit, to be investigated 

3. Allows experiments to be performed in blocks 

4. Allows designs of higher order to be built up sequentially. 

5. Provides an internal estimate of error. 

6 Does not require a laige number of runs 

7 Does not require too many levels of the independent variables 

8 Ensures simplicity of calculation of the model parameters. 

These features are sometimes conflicting, so judgement must often be applied in design 
selection 

6.4 Designs for Fitting the Second-Order Model 

A design by means of which obseiwed values of the response are collected for estimating the 
parameters in the second-order model ( Eq 6 1) is called a second-order design (Khuri and 
Cornell, 1987) 

+T.P.X +T,'LP,x.x, +s ( 6 . 1 ) 

1=1 1=1 ;=2 

where Xi, X 2 , . , Xk are the input variables which influence the response Y, Po, Pi ( i = 1,2, ..., 
k), p,j (i =1, 2, , k, j = 1, 2, .,, k ) are unknown parameters, and e is a random error 

(contribution of factors not included in (6 1)) 



Since model (6 1) contains pure quadratic terms, any second-order experimental design must 
involve at least three levels of each input variables. In the following sections various second- 
order designs are discussed 

6.4.1 The 3*^ Factorial Design 

One possible second-order design is the 3*^ factorial design, which requires that the response 
be observed at all possible combinations of the levels of k input variables which have three 
levels each In this case of 3*' factorial designs, the number N, of experimental trials is N - 
and can thus be excessively large, especially when a large number of input variables are under 
study. To reduce the total number of experimental design points, the use of fractional 
replications of these designs can be considered (Khuri and Cornell, 1987). 

DeBaun (1959) introduced a number of three-factor, three-level designs These designs consist 
of a combination of the following subsets of points taken from the 3^ design. 

1 . The 3^ factorial design 

2 Cube (points at (±1, ±1, ±1)) 

3 Center point (0,0,0). 

4 Octahedron (points at (±1, 0 , 0, 0, ±1, 0, 0, ±1)) 

5. Cuboctahedron (points at (±1, ±1, 0, ±1, 0, ±1; 0, ±1, ±1)) 

The designs investigated by DeBaun were' 

1 The 3^ factorial design 

2. Cube + octahedron + n center points, 

3 Cube + 2 octahedral + n center points. 


4 Cuboctahedron + n center points 



5. Cube + cubeoctahedion + // centei points 
6 Cuboctahcdron + octahcdion + // ccntei points 

A 

DcBaun conipaicd the above six designs on the basis ofth value {Vai(y (x))/ a^}’VN 
at (0, 0, 0) in addition to its distribution on p = (x along certain radii He concluded that 
the 3^ design is by no means the most efficient of the cases 1-6, he found that 3^ design is, in 
particular, excelled by the ci/de plus two octahedral design with 22 and 24 points, 
respectively, and by the cuboctahedron design with 16 points. 

4.4.2 Box-BehcnUin Designs 

Factorial designs for the estimation of the parameters in a second-order model was developed 
by Box and Behenkin (1960) By definition, a three level incomplete factorial design is a 
subset of the factorial combinations from a 3'" factorial design The Box-Behenkin designs are 
formed by combining two-level factorial designs with balanced incomplete block designs in a 
particular manner The interesting feature of these designs is these are “corner free”. No runs 
are done at the design corners. There are no experiments where at least one of the factors is 
not at its midpoint In contrast, to the CCD designs (discussed in the next section), there are 
no axial points so each factor appeals at only three (not five ) levels. 

6.4.3 Central Composite Designs(CCD) 

Box and Wilson (1951) introduced an alternative class of designs to the 3'^ factorial designs, 
namely, the class of central composite designs (CCD) CCD designs are considered to be most 
popular class of designs used for fitting the second-order models (Montogomery, 1997). A 
central composite design consists of 



1. A complete (or fi action of a) 2^ factoiial design, where the factor levels are coded to the 
usual -1, +1 values This is called the fadonal portion of the design. 

2 iio center points (iio ^1) 

3 Two axial points on the axis of each design variable at a distance of a from the design 
center This portion is called the axial portion of the design 

The total number of points is thus N = 2*" + 2k; + no There are two parameters in the design 
that must be specified, the distance a of the axial runs from the design center, and the number 
of center points no 

A central composite design is made rotatable by the choice of a The value of a for rotatabihty 
depends on the numbei of points in the factorial portion of the design; in fact, a - (///)*'"* yields 
a rotatable central composite design where iif is the number of points used in the factorial 
portion of the design, A choice of a in the central composite design is dictated primarily by the 
region of interest When the region of interest is sphere, the design must include center runs to 
provide reasonably stable variance of predicted response. Generally, three to five center runs 
are recommended 

6.5 Least square estimation of the regression coefficients 

The method of least squares is typically used to estimate the regression coefficients in a 
multiple linear regression model. Least square estimates minimise the sum of squares of the 
deviation between the model and the data. The method tries to minimise the sum of squares of 
the residual A model with n>k observations (k = number of independent variables {x,}, y, 
denoting the ith observed response and Xy denoting the ith observation or level of regressor Xj 


is as shown below 



yi = Po + PlX-^.l + p2X.Ki2 + 


+ Pk X ajj;. + Si 


Here Si is the error Assumptions made about the errors are(khuri) 

1 Random errors s, have zeio mean and common variance, a^. 

2 Random errors s, are mutually independent in the statistical sense 

3 Random errors s, are normally distributed. 

The method of least squares chooses the P’s in the equation so that the sum of the square of 
the errors, Si, is minimised. The above equation can be expressed in the matrix notation as 

y = Xp+6 (6.5.1) 

where 
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Xll 

••• X\k 
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• • • X2k 
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• • • Xnk 
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, and S = 
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8n 


The sum of squares of the errois has to be minimised, because the some of the error terms 
might be positive and some might be negative. The task is to minimise the deviation from the 
actual data The sum of squares of the residual can be written in the matrix form as 

z s" = (y-xp)^(y;.xp) (6.5.3) 

1 = 1 

The least squares estimate minimises the sum of squares of the residuals at the experimental 
points. Equation (6.3) is differentiated with respect to the p‘s and equating all partial 
derivatives to zero, and solved for P‘s, the least square estimates for the model is obtained. 



When the input matrix is denoted in the format given below where the vector P is the estimate 
of the parameters and x is the input variable matrix and y is the response variable vector. Then 

least square estimates, 

A 

(3 = { x^x)'^ (x^y) (6.5.4) 

The data {x, y} may be obtained from any of the statistical experimental schemes described 
earlier 

6.6 Regression models 

In this section various piocess models empirically built in the present study are presented. A 
total of four sets of second-order regression models are build. Three sets of models used the 
data generated from the MOLDFLOW PART ADVISER (PA) software and one model using 
the CMOLD 3D QUICK FILL soRware 

6.6.1 Regression models from the simulation data of PA 

As mentioned in Chapter 4 (section 4 4) operational considerations such as melt temperature 
and mold tempeiature and fill time aflect part quality by changing the conditions under which 
the mold is filled PA software allows us to change these three important process condition 
settings So, these three parameters are taken as the regressor (independent) variables and 
maximum injection pressure, temperature difference and actual injection time is taken as a 
responses. We have taken two different three factor designs for building regression models for 
the above mentioned responses. 



6.6.1. 1 Central composite rotatable design for three factors 
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Central composite rotatable design for three factors (Cochran and Cox, 1957) 


The design shown above is the central composite rotatable design for three factors, five level 
experiment This design is for estimating the parameters of the second order model. There are 
a total of twenty experiments, including six centre level experiments. Simulation experiments 
were earned out at all the twenty experimental runs Responses “injection pressure”, “actual 
injection time” and “temperature difference” were observed Each of the three responses was 
then separately used to develop the respective second-order models for injection pressure, 
actual injection time and temperature difference Least square method is used to develop the 
various response models These are shown below 

Models developed for the part “diskeady” using Polystyrene as material 
INJECTION PRESSURE = -189.01 + 0.573813 * X, ± 2.501543 * X2 -17 356 * X3 - 

0.00649 * Xi^ - 0.00695 * Xj^ + 6.184467 * X3H 0 000177 * Xi * X2 
- 0 09164 * Xi * X3 - 0 00212 * X2* X3 (6.5.1) 

ACTUAL INJECTION TIME = 0 101043 - 0 00226 * Xi - 0 00049 * X2 ± 1.06168 * X3 - 

8 2E-06 * Xi^ - 3 5E-07 * - 0.00266 * Xj^ ± 1.41E-05 *Xi * X2 


± 0 000189 * Xi * X3 ± 7.31E08 *X2 * X3 


(6.5.2) 



TEMPERATURE DIFFERENCE = 2.3741 17 + 0.03884 * X] - 0.02673 * X2 + 0.852161 * 

X 3 - 3 4E-05 * X,^ + 0 000123 * X 2 ^ - 0 5303 * Xs^ - 0 00027* X, * 
X 2 - 0 01553 * X,* X 3 + 0.01963 *X 2 *X 3 (6.5.3) 

where 

Xi = Mold temperature 
X 2 = Melt temperature 
X3 = Injection time 

Correlation diagrams for the above models relating actual values to the values obtained from 
the models are shown below 



Figure 6.6.2 
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Figure 6.6.3 


Correlation coefficients for the above models are shown in the table below 


Table 6.1 Correlation coefficients for the models developed using the design from the 
text book Cochian and Cox (1957). 


Correlation CoerTicicnts 

Injection pressure 

0 825857 

Actual injection time 

0 999995 

Temperature difference 

0 989965 


6,6. 1.2 3^ factorial design 
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3x3x3 factorial design (Box and Draper, 1987) 




The design shown above is the factorial design for a three factors, three level experiment. This 
design is for estimating the parameters of the second order model There are a total of twenty 
seven experiments in this design. Simulation runs were carried out at all the twenty seven 
experimental settings using MOLD FLOW PART ADVISER software mold the article 
“Diskcady” “Injection pressure”, “actual injection time” and “temperature difference” were 
observed as responses Each of the three responses was then separately used to develop the 
respective multiple regression models for injection pressure, actual injection time and 
temperature difference Least square method is used to develop these models for the above. 
This procedure is used to develop the second order model for the “diskcady” part using two 
different plastics, Polystyiene and LDPE. Models for both the materials are shown below. 

Models developed for the part “diskcady” using Polystyrene as material. 

INJECTION PRESSURE = 288 139 - 0 15742 * Xi - 1.45338 * Xj - 66 66 * X3 - 0.00065 

* Xi^ + 0.001644 * Xi + 7.968747 * Xs^ + 0 000945 * Xi * X2 - 
0.04881 * X, * X3 + 0 171098 * Xz^Xs (6.5.4) 

ACTUAL INJECTION TIME = 0 247042 - 0 00019 * Xi - 0 00235 * X2 + 1.06185 *X3 - 

1 8E-06 * Xi^ + 5 56E-06 * X2^ - 0 0016 * + 1.67E-06 *Xi * X2 

+ 6 68E-05 * X) * X3 + 2 22E-15*X2*X3 (6.5.5) 

TEMPERATURE DIFFERENCE = -5 87416 + 0 018823 * Xi + 0 05786 * X2 +1.363 18 * 

X3 - 5.5E-05 * Xi^ - 0 00013 * Xj^ - 0.65733 * X^^ - 8.8E-05* Xi * 

X2 - 0 01985 * Xi * X3 + 0 02146 * X2*X3 (6.5.6) 

where 


Xi = Mold temperature 
X2= Melt temperature 
X3 = Injection time 



Correlation diagrams for the above models relating actual values to the values obtained from 
the models are shown below 



Figure 6.6.5 
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Figure 6.6.6 


Correlation coefficients for the above models are shown in the table below. 


Table 6.2 Con elation coefficients for the models developed using the 3^ design for 
“Diskeady” model using Polystyrene as material 


Correlation Coefficients 

Injection pressure 

0 994314 

Actual injection time 

0 999997 

Temperature difference 

0 999188 


Models developed for the part “diskeady” using LDPE as material: 

INJECTION PRESSURE = 186.3289 - 0.14547 * Xj - 0 551 *X2- 87.2293 * X3 + 

0 000876 * Xi^ + 0 000168 * Xj^ + 12 0292 * Xa^ + 0 000244 * Xi * 
X2 “ 0 021 13 * X, * X3 + 0. 149972 * X2 * X3 (6.5.7) 

ACTUAL INJECTION TIME = -0 21424 + 1 36E-05 * Xi + 0.001889 * X2 + 073517 * X3 

- l.OE-16 * Xj^ - 4 2E-06 * Xa^ + 0 000746 * Xj^ -3 3E-16 * Xi * 


X2 + 2 79E-5 * Xi ^ X3 - 2 8E-05* X2 * X3 


(6.5.8) 





TEMPERATURE DIFFERENCE = -4 75201 - 0 00028 * Xi + 0 045121 * Xz + 2.40199 * 


Xs - 4 OE-05 * X,^ - 9 OE-05 * Xz^ - 0.79691* - 1.7E-05* Xi * 

Xz - 0 02026 * Xi * X3 + 0.01847 * Xz * X3 ( 6 . 5 . 9 ) 

where 

Xi = Mold temperature 
Xz = Melt temperature 
X 3 = Injection time 

Correlation diagrams for the above models relating actual values to the values obtained from 
the models are shown below 



Figure 6.6.8 
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Figure 6.6.9 


Correlation coefficients calculated for the above models are shown in the table below. 
Table 6.3 Correlation coefficients for the models developed using the 3^ design for the 
“Diskcady” model using LDPE as material. 


Correlation Coel 

fficients 

Injection pressure 

0.997018 

Actual injection time 

0.999997 

Temperature difleience 

0.999137 


6.6.2 Regression models from the simulation data of 3D Quick Fill 

This software allows to change the process conditions melt and mold temperature and the 
molding machine parameters used for the process, the maximum injection pressure and the 
maximum injection rate. Though the machine parameters are not real process conditions, we 
have taken these parameters also as regressor (independent) variables as the 3D Quick Fill so 
showing change in the responses by changing these parameters Injection pressure, fill time, 
cooling time and temperatuie difference are taken as responses As there are four independent 
variables we have used Box and Behenkin four factor design for building the regression 


models 












6.6.2. 1 Box and Behenkin four factor design 

'±1 ±1 0 o' 

0 0 ±1 ±1 
0 ±1 0 0 
±1 0 ±1 0 
0 0 0 ±1 
0 0 0 0 
0 ±1 ±1 0 
±10 0 ±1 
0 0+10 
±10 0 0 
_ 0 ±1 0 ± 1 _ 

Box and Behenkin Design for Four Factors (Jugal, 1997) 

Shown above is the Box and Behenkin design for a four factor, three level experiment. Total 
thirty three runs are performed using CMOLD 3D QUICK FILL software for “tape” model 
Parameters are estimated from the generated data for second order response model. Least 
square parameter estimate method is used to estimate the parameters Responses for this 
experiments are “injection pressure”, “fill time”, “cooling time” and “temperature difference”. 
Regression models are developed for all four responses are listed below. 

INJECTION PRESSURE = 125.494 - 0.590052 * Xi - 0.03866 * X 2 + 0.370962 * X 3 + 

0 004097 * X 4 + 0.000786 * Xi^ - 0.00017 * X 2 - 0.00037 * - 2.1E-07 

* X,^ + 0 000206 * Xi* X 2 - 0 00073 * Xi * X 3 + 1 38E-06 *Xi*X-7E- 

05 * X 2 * X 3 - 6E-07 * X 2 * Xj + 6.07E-06*X3 * X 4 , (6.5.10) 

FILL TIME = 0 013339 - 2. lE-07 * Xi - 2 5E-07 * X 2 - 2.4E-05 ♦ X 3 - 1.3E-06 * X + 

1 99E-08 * Xi^ - 8 9E-09 * X 2 ^ + 2 34E-08 * Xa^ + 5.56E-1 1 * X^ + 

1 87E-09 * X* X 2 + 2 99E-08 * X * X 3 + 2.81E-10 * XI * X + 2 25E-09 

* X 2 * X 3 + 2 25E-1 1*X2 * X 4 + 3.89E-10*X3 * X (6.5.11) 



COOLING TIME = - 0 10767 + 0 02642 * Xi - 0 09467 * Xj - 0 00023 * Xj - 1 4E-05 * 
X4 - 3.5E-05 * X,^ + 0.001353 * X-z + 1.94E- 07 * Xs^ - 2. lE-10 * X4^ + 

1 17E-05 * Xi* X2 - 9 4E-08 * Xj * X3 + 2 81E-08 * X, * X4 - 1 88E-06 * 
X2 * X3 - 3 75E - 08* X2 * X4 + 7.5E-09 * X3 * X4 ( 6 . 5 . 12 ) 

TEMPERATURE DIFFERENCE = 8 899248 - 0 04484 * Xi + 0 009147 * X2 - 0.00991 * 
X3 - 0.00034 * X4 + 5 89E-06 * Xi^ - 7E-05 * Xz + 5.19E-06 * Xs^ + 

9 17E-09 * X4^ - 4 lE-15 * Xi* Xi + 2 81E-05 * Xi * X2 + 5.62E-07 * XI * 
X4 - 4.3E-17 * X2 * X3 - 3,2E-13 * X2* X4 - 1.3E-13 * X3 * X4. ( 6 . 5 . 13 ) 

Where 

Xi. Melt temperature 
X2 Mold temperature 
X3: Maximum injection pressure 

X4. Maximum injection rate 

Correlation diagrams for the above models relating actual values to the values obtained from 
the models are shown below 
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Figure 6.6.13 

Correlation coefficients are calculated for the above models, which are shown in the table 
below 


Table 6.4 Correlation coefficients for the models developed using the Box and Behenkin 
design 


Correlation Coefficients 

Injection pressuie 

0.997001 

Fill time 

0 999460 

Cooling time 

0.999995 

Temperature difference 

0 977995 


Summary 


This chapter has developed the multiple 'regiession process models that are subsequently used 
to optimise the responses by GA. 




Chapter 7 


RESULTS FROM SGA AND NSGA 

In this chapter results obtained from the different single and multi-criteria optimisation 
problems attempted in this study are presented The simple genetic algorithm (SGA) and non- 
dominated sorting genetic algorithm (NSGA), a new meta-heuristic are used for optimisation. 
Both SGA and NSGA are applied to four sets of empirically second-order regression models 
developed. Convergence graphs are shown from the results obtained from the execution of 
SGA. Paieto diagrams resulting from multi-criteria optimisation by NSGA are shown. Solution 
populations after fifty NSGA generations for several different multi-criteria optimisation 
attempted are shown in tables 

As mentioned in the Section 5 2 2, the Moldflow simulator Part Adviser allows one to change 
three process conditions, namely, “melt temperatuie”, “mold temperature” and “injection time”. 
So, we used these three variables as process variables for building the regression models From 
the various results that are given by Part adviser (listed in Section 5.2.2) we took “injection 
pressure’, “temperatuie difleience” (noted from flow front temperature result) and “actual 
injection time” as i espouses. 

The CMOLD simulator 3d Quickfill allows one to specify the process conditions “mold 
temperature” and “melt tempeiature” It also allows one to specify three other molding machine 
parameters for the process, namely, the maximum injection pressure, the maximum injection 
rate and the machine performance. So, we used mold temperature, melt temperature, maximum 
injection rate and maximum injection pressure as process variables for developing regression 



models For the various results that 3d Quickfill (listed in the Section 5.1.2) given we took 

“injection pressure”, “fill time” , “cooling time” and “temperature difference” as responses. 

\ 

7.1 Evaluation of Factor effects of Moldflow Simulator responses 

To get the insight into which factors are important in a particular process one has to do the 
screening (David and Robert, 1998) One would accomplish this is by running a given factor at 
two levels (a high level and a low level) within a DOE framework and seeing if varying the 
level of this factor had any effect on the response. The simplest experimental design for 
accomplishing this is the Two level full factorial design So, to get the insight into responses of 
the Moldflow simulatoi, experiments were run using a 2^ full factorial design The effect of 
varying the three process variables of Moldflow simulator on the three responses are shown in 
the Figures 7 1 to 7 3 

Engineering considerations tell that increasing the melt and mold temperatures will reduce the 
required injection pressure, which we can see in Figure 7.1. As the melt and mold 
temperatures increase, the melt viscosity decreases hence there will be decrease in the injection 
pressure However, in the relationship between injection time and injection pressure, several 
contradictory factors come into play The shorter the injection time the higher is the 
volumetric flow rate, and the higher the pressure requirement. However, high speed injection 
also generates frictional (viscous heating) that raises the material temperature. The combined 
effect of high temperature and high shear rate (resulting from high flow rate) reduces the melt 
viscosity, and therefore offsets the pressure requirement. In the experiments see from the 
Figure 7. 1 that with increase in injection time there is slight increase in injection pressure. 



Figure 7.2 Factor effects on Temperature dilTerence, Moldflow Simulator, 2^ DOE 



Figure 7.3 Factor effects on Actual injection time, Moldflow Simulator, 2^ DOE 



Mold Temp 


Melt Temp 


Inj.Time 







7.2 Results from Executing the SGA 


When we use the SGA, we are interested in optimising a single objectives and in finding the 
single best solution GA responses chosen are usually the best, average, and the worst value of 
the objective function found after executing a fixed number of GA iterations Second-order 
regression models developed for each of the responses may be separately plugged in SGA for 
response evaluation These included the “injection pressure” , “actual injection time” and 
“temperature difference” regression models built using the simulation runs of Moldflow 
simulator “Injection pressure”, “fill time”, “cooling time” and “temperature difference” 
regression models built using the simulation runs of CMOLD simulator were used for the 
CMOLD processes 

A fundamental result in genetic algorithms, the schema theorem (Goldberg, 1989), by itself is 
insufficient to guide the effective use of GAs in global optimisation It does not indicate what 
specific values should one choose for the different GA parameters, namely Population size. Pc 
and Pm etc , and the different schemes that introduced controlled amount of randomness while 
the GA is executing(Jugal, 1997) However, it is well known now that the GA’s efficiency 
depends on high degree upon the selection of these “control parameters’ (Davis, 1991) A 
further complication exists because it is also widely reported that “optimum parameter settings” 
may be problem-specific, and that the effects of these parameters may interact Indeed, Davis 
notes that crossover and mutation effects can interact and “support each other in important 
ways” (Davis, 1991) Davis observes that a judicious blend of mutation and crossover does 
better than either one alone to strike a good balance between exploration of the total solution 
space and exploitation of good solutions currently at hand As it is noted above, the successful 
implementation of GA would depend upon the correctness in the choice of different parameters 



that control its execution A two level design of experiments is conducted to find the best 
combination of Pc and Pm keeping the other GA parameters constant. We found that the 
combination of 0.9 crossover probability and 0.005 mutation probability is giving the best 
results for hundred population So following combination of GA parameters are used for SGA 
application' 

Population size = 100 

Probability of crossover = 0 9 

Probability of mutation = 0 005 

Table 7.2. 1 to Table 7.2 6 shows the best, average and worst values of responses after 50 
generations of a four different experiment runs for Process models developed using 3^ design 
on moldflow simulator. Table 7 2.1 to Table 7 2.3 shows the results for the second-order 
regression models for polystyrene diskcady and Table 7.2.4 to Table 7.2.6 is for LDPE 
diskcady. 

Table 111 Xo Table 7.2.9 shows the best, average and worst values of responses after 50 
generations for Process models developed using the CCD rotatable design on moldflow 
simulator for polystyrene diskcady 

Table 7 2.10 to Table 7.2.13 shows the best, average and worst values of responses after 50 
generations for Process models developed using the Box and Behenkin design on CMOLD 
simulator for polypropylene tapecasing. 

Figure 7 2 1 to Figure 7 2 26 shows the convergence graphs for best and average values for 
the various second-order regression models developed after 50 generations 



Table 7.2.1 : Injection pressure (MPa) values after 50 Generations for 
process model by 3^ DOE on Moldflow Simulator for 
Polystyrene Diskcady 


Seed 

Best 


Worst 

1 


15 6030 

70.3796 


13 7858 

14 8402 

74 2196 

EHHI 

14 0934 

15 5381 

70 3709 


13 9178 

16 4729 

67 6516 


Table 7.2.2 : Actual injection lime (Sec) values after 50 Generations for 
|)roccss model by 3^ DOE on Moldflow Simulator for 
Polystyrene Diskcady 


Seed 

Best 

iroiWHWII 


1 

0 0102 




0.0128 

0 0708 

3.0334 

50 


0 0489 

3 0326 

75 


0 0509 

3 0681 


Table 7.2.3 : Temperature difference (°C) values after 50 Generations for 
process by 3^ DOE on Moldflow Simulator for Polystyrene 
Diskcady 


Seed 

Best 

Average 

Worst 

1 

0 2826 

0 3542 

10 224 

30 

0 2826 

0 3875 

10 224 

50 

0.2901 

0.4042 

10 224 

75 

0 2826 

0 4056 

10 224 

















Table 7.2.4 : Injection pressure (MPa) values after 50 Generations for 

process model by 3^ DOE on Moldflow Simulator for LDPE 
Diskcady 


1 Seed 

Best 

Average 

Worst 

WjHIHII 

3 7623 

4 2558 

69 5169 

30 


4.0710 

69 5169 

50 

3 7529 



75 



69 5169 


Table 7.2.5 : Actual injection time (Sec) values after 50 Generations for 
process model by 3^ DOE on Moldflow Simulator for LDPE 
Diskcady 




Average 






1 30 

001 

0 0772 

3.0619 

jEEHIlii 

001 

0 0588 

3.0585 


0 01 

0.0673 

3.0924 


Table 7.2.6 : Temperature difference (°C) values after 50 Generations for 
process by 3^ DOE on Moldflow Simulator for polystyrene 
Diskcady 


Seed 

Best 

Average 

Worst 

1 

0 361 

0.4482 

11.568 

30 

0.361 

0.4501 

11.568 

50 

0 361 

0 4511 

11 568 

75 

0 361 

0 4193 

11.568 

















Table 7.2.7 : Injection |)ressiire values (MPa) after 50 Generations for 
process model by CCD rotatable design on Moldflow 
Simulator for Polystyrene Diskcady 



Best 

Average 

Worst 

1 

9 6293 

11 091 

46 8810 

30 

94514 

10.4182 

49.334 

50 

9 4167 

10 3166 

45 7444 

75 

9 6317 

11.7011 

46.9483 


Table 7.2.8 : Actual injection time (Sec) values after 50 Generations for 

process model by CCD rotatable design on Moldflow Simulator 
for Polystyrene Diskcady 


1 Seed 

Best 

Average 

Esmmm 

I^^Hi 





0 0075 

0 431 


1 50 

0 0075 

0 532 

3.0314 



0.0433 

3.0586 


Table 7.2.9 : Temperature difference (°C) values after 50 Generations for 
process model by CCD rotatable design on Moldflow 
' Simulator for Polystyrene Diskcady 


Seed 

Best 

Average 

Worst 

1 



10 2395 

30 

0 7446 


10.8202 

50 

0 7466 


10.2068 

75 

0 7446 

0 8329 

10 2068 




























Table 7.2.10 : liijectiun pressure values (MPa) after 50 Generations for 
process model by Box and Belienkin design on CMOLD 
Simulator for Polypropylene Tapecasing 


Seed 

Best 

Average 

Worst 

1 

33 2792‘ 

37 1807 

105 9656 

30 

32 8074 

36 7974 

105 9656 

50 

32 665 

36 4936 

105 9656 

75 

32 4539 

36 6764 

105 9656 


Table 7.2.11 : Fill time values (Sec) after 50 Generations for process model 
by Box and Belienkin design on ^MOLD Simulator for 
Polypropylene Tapecasing 


Seed 


Average 

Worst 

1 

0 0010 


0 0085 





50 

0.0010 



75 

0 0010 

0 0026 

0 0085 1 


Table 7.2.12 : Cooling time values (Sec) after 50 Generations for process 
model by Box and Behenldn design on CMOLD Simulator 
for Polypropylene Tapecasing 


Seed 

Best 

Average 

^ESSSQHHI 


2.1566 


BBESHI 

30 

2.1560 


5 7449 

50 

2 1539 

2 4354 

5 8153 

75 

2 1570 

2 3984 

5 7449 


Table 7.2.13: Temperature difference (°C) values after 50 Generations for 
process model by Box and Belienkin design on CMOLD 


S imulator for Polypropylene Tapecasing 


Seed 

Best 

gw 

Worst 

1 

0 1527 

EESdHi 

2 0417 

30 

0 1490 

BS9BB 

1.9269 

50 

0.1561 


1 9269 

75 

0 1506 

0 1933 

1 9629 
































Figure 7.2.1 Convergence graph 

3x3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.2 Convergence graph 

3x3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.4 Convergence graph 

3x3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 









Figure 7.2.5 Convergence graph 

3x3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.6 Convergence graph 

3x3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.7 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Figure 7.2.8 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 





Figure 7.2.9 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Figure 7.2.10 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 





Figure 7.2.11 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskc 
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Figure 7.2.12 Convergence graph 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 
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Figure 7.2,13 Convergence graph 
CCD model, Moldllow Simulator, Polystyrene Diskcady 
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Figure 7.2.14 Convergence graph 
CCD model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.15 Convergence graph 
CCD model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.16 Convergence graph 
CCD model, Moldflow Simulator, Polystyrene Diskcady 


e 

Q 

s 

I 

I 

H 

H 

CO 

pa 


u 


(5 

5 

4 

3 

2 

1 

0 l_ 

1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 

GENERATIONS 



Figure 7.2.17 Convergence graph 
CCD model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7;2.18 Convergence graph 
CCD model, Moldflow Simulator, Polystyrene Diskcady 
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Figure 7.2.19 Convergence graph 

Box-Behenkin model, Moldllow Simulator, Polypropylene Tapecasing 



Figure 7.2.20 Convergence graph 

Box-Behenkin model, Moldllow Simulator, Polypropylene Tapecasing 





Figure 7.2.21 Convergence graph 

Box-Behenkin model, Moldflow Simulator, Polypropylene Tapecasing 



Figure 7.2.22 Convergence graph 

Box-Behenkin model, Moldflow Simulator, Polypropylene Tapecasing 





Figure 7.2.23 Convergence graph 

Box-Behenkin model, Moldflow Simulator, Polypropylene Tapecasing 



Figure 7.2.24 Convergence graph 

Box-Bchcnkin model, Moldflow Simulator, Polypropylene Tapecasing 




Figure 7.2.25 Convergence graph 

Box-Behenkin model, Moldflow Simulator, Polypropylene Tapecasing 
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Figure 7.2.26 Convergence graph 

Box-Behenkin model, Moldflow Simulator, Polypropylene Tapecasing 
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7.3 NSGA results 


In this study total three experimental designs used to build four sets of regression models 
Three from the data of Moldflow simulator and one fiom the data from the CMOLD simulator 
Various designs build and their details are as below 

1 3^ full factorial DOE is used to develop the second-order regression models for the three 
responses of the Moldflow simulator for the Polystyrene Diskcady and LDPE Diskcady 

2 CCD rotatable design is used to develop the second-order regression models for the three 
responses of the Moldflow simulator foi the Polystyrene Diskcady 

3 Box and Behenkin four factor design is used to develop the second-order regression models 
for the foul responses of the CMOLD simulator for the Polypropylene Tapecasing 

Various second-order regiession models developed are shown in the section 6,6 

There are totally four sets of regression models, three sets containing three second-order 
regression models each (Moldflow simulator) and one set contains four second-order 
regression models (CMOLD simulator) Multi-criteria optimisation is attempted for these 
second-order regiession models using NSGA We have considered two objectives at a time for 
simultaneous mminiisalion As iheie aie ihice losponses foi the Moldflow simulator total ol 
three Bi-objective problems for each set of regression models The three Bi-objective problems 
are 

• Injection pressure and Actual injection time 

• Injection pressure and Temperature difference 

• Temperature difference and actual injection time 

As there are four responses for the CMOLD simulator total six bi-objective problems are 
tested Vaiious paicto diagiams resulting fiom the optimisation of the above problems arc 



shown and the pareto solutions after fifty generations are shown in the tables The six bi- 
objective problems are 

• Injection pressure and Fill time 

• Injection pressure and Cooling time 

• Injection pressure and Temperature difference 

• Fill time and Temperature difference 

• Fill time and Cooling time 

• Temperature difference and Cooling time 

7.3.1 Parameters used in NSGA Execution 


Following GA parameter values and the limits of process variables used in execution of 
NSGA for the Bi-objective problems of Moldflow simulator for Polystyrene Diskcady A two 
level DOE is conducted to find the best combination of Pc and Pm keeping the other NSGA 
parameters constant We found that the combination of 0 9 crossover probability and 0 01 
mutation probability is giving the best results for hundred populaion 


• Population size = 100 

• Per OSS (Pc) = 0 9 

• Pmute (Per) ~ 0 01 

• Dshare = 1 

• Length of chromosome = 24 

• Mold Temperature (mm) = 20 °C 

• Mold Temperature (max) = 70 °C 

• Melt Temperature (mm) = 185 °C 

• Melt Temperature (max) = 225 °C 

• Injection time (min) = 0 0075 Sec 

• Injection time (max) = 3 Sec 



Following GA parameter values used and the limits of process variables in execution of NSGA 
for the Bi-objective problems of Moldflow simulator for LDPE Diskcady 


0 

Population size 

100 


0 

Pcross (PJ 

09 


0 

Pmute (P„J 

001 


0 

Dshare 

1 


0 

Length of chromosome = 

24 


0 

Mold Tempeiatw e (min) = 

20 °C 


0 

Mold Tempetature (max) = 

60 °C 


0 

Melt Tempeiatiiie (min) - 

205 °C 

0 

Melt Tempeiatiire (max) = 

245 °C 

0 

Injection time (min) = 

0 01 Sec 

0 

Injection time (max) = 

3 Sec 


Following GA parameter values and the limits of process variables that are used in NSGA 

execution for the Bi-objective problems of CMOLD simulator for Polypropylene Tapecasing 

0 

Population size 


100 


Pcross (Pf 


09 


Pmute (Pm) 


0 01 


Dshare 


1 


Length of chromosome 

= 

32 


Mold Temperature (min) 

= 

20 °C 


Mold Temperature (max) 


00 

o 

o 

n 


Melt Tempetature (min) 


200 °C 


Melt Tempetature (max) 


280 °C 


Max Injection Pressure (mm) 


0 Mpa 


Max Injection Pressure (max) 

= 

300 Mpa 


Max Injection Rate (min) 

=- 

0 cm^/Sec 


Max. Injection Rate (Max) 

- 

10000 cm^/Sec 
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Figures 7 3 1 show the initial population of chromosomes and the population after fifty 
generations Figure 7 3 2 show the movement of the members as NSGA executes on to the 
non-dominated (Pareto) fi’ont for every fifth generation from randomly generated initial 
population to twenty fifth generation Both these figures is for the Bi-objective problem of 
“injection pressure” and “actual injection time” for polystyrene diskcady Table 7 3 1 show the 
Pareto optimal solutions after fifty generations Subsequently from the resulting solutions 
(Table 7 3 1) on the Pareto front the decision maker will be expected to select any solution 
depending his preference criteria and performance requirements A particular specification of 
Mold temperature, Melt temperature and Injection time yields the corresponding values of the 
Injection pressure and actual injection time Figures 7 3 3 to Figures 7 3 20 show the Pareto 
solutions obtained for the other Bi-objective problems 

Figures 7 3 13 and Figure 7 3 14 show the solutions obtained for the Bi-objective mimmisation 
problem of “actual injection time” and “temperature difference” after 50 generations In this 
figures we can notice that there is no Pareto fi'ont and the solutions converged at the low left 
comer This can be explained as follows with increase in injection time there is increase in 
temperature difference along the flow fi'ont As objectives or not conflicting there is no 
possibility of Pareto fi'ont. As we are doing minimisation of the two objectives simultaneously 
the solutions converged to the low left comer of the pareto diagram The same observation 
can also be noted in the Figure 7 3 17 for the Bi-objective problem of “fill time” and 
“temperature difference”, which is from the results of CMOLD simulator 
Tables 7 3 2 to Tables 7 3 14 show the Pareto optimal solutions that are obtained from the 
execution of NSGA for the various Bi-objective problems 

The results that are obtained from the NSGA are cross checked by using the process variable 
values given by the NSGA for the different solutions along the pareto front for all the cases 



Actual injection time (Sec) Actual injection time (Sec) 


Figure 7.3.1 Pareto solutions after 50 generations 
3x3x3 DOE model, Moldflow simulator, Polystyrene Diskcady 



Injection pressure (MPa) 


Figure 7.3.2 Convergence of solutions to Pareto front 
X 3 X 3 DOE model, Moldflow simulator, Polystyrene Diskcady 



Injection pressure (MPa) 
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Injection pressure (MPa) w Injection pressure (MPa) 



Figure 7.3.4 Convergence of solutions to Pareto front 
3x3 DOE model, Moldflow Simulator, Polystyrene Diskcady 
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Actual injection time (Sec) Actual injection time (Sec) 


Figure 7.3.5 Pareto solutions after 50 generations 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Injection pressure (MPa) 


Figure 7.3.6 Convergence of solutions to Pareto front 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Injection pressure (MPa) 
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Injection pressure (MPa) 


Figure 7.3.7 Pareto solutions after 50 generations 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Temperature difference (^) 


Figure 7.3.8 Convergence of solutions to Pareto front 
3x3x3 DOE model, Moldflow Simulator, LDPE Diskcady 



Temperature difference (^) 
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Actual injection time (Sec) Actual injection time (Sec) 


Figure 7.3.9 Pareto solutions after 50 generations 



Injection pressure (MPa) 


Figure 7.3.10 Convergence of solutions to pareto front 
CCD model, Moldflow Simulator, Polypropylene Diskcady 
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Injection pressure (MPa) 


Figure 7.3.11 Pareto solutions after 50 generations 



Temperature difference ( C) 


Figure 7.3.12 Convergence of solutions to Pareto front 


CCD model, Moldflow Simulator, Polypropylene Diskcady 



Temperature difference ( C) 
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Actual injection timeJSec) ^ Actual injection time (Sec) 


Figure 7.3,13 Pareto solutions after 50 generations 
3x3x3 DOE model, Moldflow simulator, Polystyrene Diskcady 
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Figure 7.3.14 Convergence of solutions to Pareto front 
X 3 X 3 DOE model, Moldflow simulator, LDPE Diskcady 



Temperature difference (®C) 



132 




Figure 7.3.15 Pareto solutions after 50 generations 
X and Behenkin design model, CMOLD Simulator, Polypropylene Tapecasing 



Figure 7.3.16 Convergence of solutions to Pareto front 
c and Behenkin design model, CMOLD Simulator, Polypropylene Tapecasing 
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Figure 7.3.17 Pareto solutions after 50 generations 
Box and Behenkin model, CMOLD Simulator, Polypropylene Tapecasing 



Temperature difference ( (?) 


Figure 7.3.18 Convergence of solutions to Pareto front 
Box and Behenkin model, CMOLD Simulator, Polypropylene Tapecasing 
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Figure 7.3.19 Pareto solutions after 50 generations 



Injection pressure (MPa) 


Figure 7.3.20 Convergence of solutions to Pareto front 
Box and Behenkin design model, CMOLD simulator, Polypropylene Tapecasing 
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Following notations are used in the Tables 


Mold Temp 
Melt Temp 
Inj Time 
Inj Pres 
Ac Inj Time 
Temp DifF 
M Inj Pre 
M Inj. Rate 


= Mold Temperature 
= Melt Temperature 
= Injection time 
= Injection Pressure 
= Actual Injection Time 
= Temperature Difference 
= Maximum Injection Pressure 
= Maximum Injection Rate 


Table 7.3.1 : NSGA population for process model by 3^ factorial design 


Moldflow Simulator for Polystyrene Diskcady 


[ Pareto Ootimal Solutions after 50 Generations 1 

I Process variables 

Responses 

Mold Temp. 

°c 

Melt Temp. 

"C 

Inj. Time 

Sec 

Inj. Pres. 

MPa 

Ac.Inj.Time 

Sec 

69 41 

223.75 

1 52 

16 01 

1 62246 

69 41 

222 18 

1 05 

22 04 

1 123387 

68 63 

224 84 

1 86 

14 11 

1 983806 

68 63 

202 1 

0 04 

5932 

0 046388 

66 27 

224 84 

1 87 

14 38 

1 996091 

69 8 

225 

1 54 

15.31 

1.647718 

69 02 

214.96 

1.54 

19.41 

1 645565 

69 41 

224 84 

1.71 

14 48 

1 822024 

69 8 

224 84 

1 83 

14 03 

1 946537 

69 41 

224.69 

Oil 

41 8 

0.124522 

57 25 

225 

1 17 

20 1 

1.248401 

69 8 

219 98 

1.64 

16 61 

1 746159 

69 41 

224 84 

1 87 

13.99 

1 996299 

66 27 

225 

2 08 

14 35 

2 22003 

69 8 

224.69 

0 66 

27 76 

0.712058 

69 8 

225 

1 64 

14 73 

1 747367 

698 

225 

1 51 

15 56 

1 610343 

698 

225 

02 

39 2 

0 212154 

43 53 

204 92 

1 93 

24 3 

2053405 

69 8 

225 

1.15 

19 35 

1.223897 

69 8 

224 84 

1.83 

14 03 

1.946537 

69 8 

224.69 

066 

27.76 

0.712058 
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69 8 


68 63 


69 41 


70 


69 41 


69 8 


69 8 


69 41 


69 41 


69 8 


69 8 


63 14 


63 53 


69 41 


69 8 


6 


69 41 


62 55 


70 


68 63 


69 41 


68 63 


68 82 


69 8 


68 24 


44 71 


69 8 


223 43 


225 


225 


224 69 


224 69 


225 


225 


224 84 


224 69 


225 


225 


224 53 


223 43 


219 98 


203 5 


31 87 


39 59 


14 08 


18 93 


224 84 


225 


224 69 


224 69 


69 41 

224 69 

66 67 

225 

31 76 

225 

69 41 

224 69 


22 13 


13 88 


43 84 


45 03 


41 94 


31 89 


22 69 


28 01 


22 01 


45 25 


43 37 


189 


44 82 


43 78 


18 45 


20 96 


0 549315 


0 199705 


1.946515 


1 273663 


0 137095 


0 724549 


1 323575 


1 024276 


1 996365 


0 049488 


0 108956 


0 011877 


0 1 12072 


0 911917 


1 273505 


0.549328 


1 147845 


0 896318 


0.724603 


0 721489 


1 995798 


1.996334 


1.298705 


1 523023 


1.298705 


1 398367 


0 023268 


1 211337 


1 049217 


0011831 


0.674592 


0 06202 


1 273699 


0 011961 


0.0496 


1 298726 


1 247393 


1.710002 
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44 38546 


35 43248 


63 73 


70 


57 45 


44 9 
70 


70 


44 9 


66 86 


70 


63 53 


70 


69 8 


68 43 


70 


68.43 


70 


69 8 


70 


224 84 


53 


BEoHI 



224 53 


224 84 


224 84 


222 33 


224 53 


225 


224.84 


224 84 


224.84 


223 9 


225 


22 


224 84 


224.84 


70 

224 84 

69 8 

224.84 


23 17444 


23 06239 


43 65834 


20 75549 


41 85003 


4 01 

23 17444 

4 14 

22 57623 

3 64 

28 97431 

2 52 

31 69531 

4 14 

22 57623 

3 85 i 

24 01798 

4 15 

^'92936 


33 03829 


22 20458 


22 54996 


21.86379 


30 78458 


21 75392 


30 78458 


22.85628 
26 61609 ' 


21 00357 


35 43248 


23 17444 


21 07755 


23 


21 62309 


21.63962 


24 01798 


22 97284 


224 53 


224.84 


22 


41 53798 


31 79445 
18 45788 


25 76486 


30 78458 


43 65834 
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Table 7.3.3 : NSGA population for process model by 3^ factorial design on 
Moldflow Simulator for Polystyrene Diskcady 


Pareto Optimal Solutions after 50 Generations 


Responses 


Temp. DifT. Ac.Inj. Time 





0 409253 


0 01084 


0 009351 


0 009014 


0 008798 


0 009834 


0 010478 



0 009884 


0 008817 


1 709221 


0 009096 


0.009989 


0 008952 


0 209135 


0 008757 


0 409591 


0,010042 


0 01084 


0 010772 


0 008813 


0.310943 


0 059355 



0 009936 


0 008857 


0 061074 


0 008952 


0 00916 


0 109971 




























































































70 

70 

69 22 
63 73 


57 45 


196 61 
216 22 
222 49 
"215 43 
215 12 


6 


206 65 


96 6 


216 22 


215.12 


207 59 


0 11 

0 79 

0.2 

] 

27 


0 011362 


0 122612 


0 008796 


0 010015 


0 009834 


0 008952 


70 


68 43 


7 


70 


70 


7 


7 


66 86 


69 8 


70 
69 8 
70 


8 


211 67 


215 43 


216 22 


208.84 
207 59 


220.45 


207 59 


5 59 


211 82 


95 35 


22 


2 


2 


206 18 
205 . 08 " 
220 45 


215.75 


218 57 


0 010042 


0 008798 


0.009351 


0 008952 


031 

0 009884 

031 

0 010015 


0 008952 


0 01084 


0 008952 


0 00991 


0 010854 


0011214 


0 009371 


0 008863 


0 01084 


0 010276 


0 010031 


0 34 

0 008857 

0 34 

0 008817 

04 

0 035871 

031 

0 0099'89 

0 33 

0 009055 


0 209491 


0 009936 


0 010449 


70 

207 59 

0 01 

0 34 

0 008952 

70 

205.08 

0.01 

0.34 

0 008798 
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2 


22 


205 08 


69 8 

216 22 

70 

207 59 

70 

187.51 

38 63 " 

209 31 

70 

208 84 

70 

215 75 

70 

211 35 

70 

204 14 


■ 


■ 

0 011214 

[ 

0 38 

0 009348 

1 

031 

0 010031 

[ 

0 34 

0.008952 


0,009682 


0 010211 


0 009055 


0 009936 


66 86 


70 


70 


70 


70 


208.22 


210 88 


216 22 


217 47 


5 


220 6 



0 34 

0 00894 


0 33 

0 009351 


0 009884 


0 009041 


0 01 

0 36 

0 009133 

001 

0 34 

0 009001 

001 

0 33 

0 00926 

001 


0 010015 

001 

03 

0 010239 

0.01 

0.29 

0 01084 


0 060298 


0.010875 


0 008952 


0 009936 
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Table 7.3.4 : NSGA population for process model by 3^ factorial design on 
Moldflow Simulator for LDPE Diskcady 


1 Pareto Ootimal Solutions 1 

1 Process Variables 

Responses 

Mold Teni|). 

°c 

Melt Temp. 

°C 

Inj. Time 

Sec 

Inj. Pre. 

MPa 

Ac, Inj. Time 

S/'C 

58 9 

245 

0 1 

54 17 

0 1 OS 163 

39 14 

244 84 

0 85 

24 4 

0 90'-'676 

29 57 

244 84 

1 46 

10 64 

1 561245 

33 8 

244 53 

061 

33 01 

0 646383 

38 35 

244 84 

0 24 

47 66 

0 258123 

41 18 

245 

1 49 

9 6 

1.586937 

38 98 

244 84 

0 08 

55 43 

0082818 

41 49 

245 

0 85 

24 27 

0 909734 

56 86 

245 

0 01 

58 86 

0 007916 

28 31 

224 76 

146 

15 66 

1 563674 

51 69 

244 69 

0 16 

51 4 

0 170729 

38 04 

238 73 

1 04 

20 96 

1 11126 

41 49 

244 84 

0 24 

47 6 

0 258187 

38.98 

244 69 

2 09 

4,77 

2 226992 

37 88 

244 69 

0 16 

51 55 

0,170479 

37 88 

243 59 

1 04 

19 48 

1 110447 

38.98 

244 84 

0 08 

55.43 

0 082818 

30 51 

243 59 

2 07 

5 53 

2 214072 

58 9 

245 

0,1 

54 17 

0 108163 

41 96 

239 82 

0,55 

36 72 

0 584733 

20.94 

241 71 

1 89 

7 34 

2 012984 

37 88 

244 69 

0 16 

51,55 

0 170479 

21 1 

244 84 

0 24 

48 33 

0,257771 

38 35 

245 

0 57 

33 98 

0 608846 

31 45 

245 

0 1 

54 41 

0 10771 

31 45 

245 

0 85 

24,67 

0 909358 

51 06 

245 

0,85 

24,06 

0 910092 

31 45 

244 22 

0 87 

24 58 

0,922039 

21 1 

244 84 

1 89 

6 76 

2 012344 

51.22 

244 84 

0 01 

58 88 

0 007864 

49 02 

243,43 

24 

5 12 

2.567186 

36 16 

244 69 

0 54 

35.51 

0 571252 

38 04 

244 84 

06 

33.18 

0 633925 

40 39 

245 

0,12 

53 62 

0 12038 
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41 18 
41 18 


32.08 


40 39 


41 18 


41 65 


41 18 


41 49 


52 63 


47 61 


58 9 
49 18 


37 88 


41 65 


245 
244 84 


244 84 


245 


245 


244 69 


245 


244.69 


244 53 


41,33 

244 69 

21 1 

245 

41 18 

245 

37 1 

244.22 

21 1 

244 84 


28 31 

239 98 

48.86 

244 69 

38,98 

245 

5671 

244 84 

31 14 

244 84 


33 9 
5 33 


33 4 


24 31 


34 2 


24 17 


53 73 


24 36 


58 86 


33 4 


33 07 


34 65 


34.02 


34 77 


33 9 


59 56 


0 608929 


2 013674 


0 633745 


0 909692 


0 608669 


0 922424 


40 

ri 


0 120454 


1 586937 


0 90964 


0 007702 


245 

0 57 

34 01 

0 608818 

245 

1 09 

18 

1 160183 

244.84 

09 

22 92 

0 95989 


0 62186 


0 634272 


0 973125 


0 5841 


0.608991 


0 608335 


0 007446 


51 06 


41 18 


35 3 

0 571363 

22 9 

0 959908 

33 48 

0 634236 

18 78 

1 110621 

56 77 

0 108415 

35 67 

0 559083 

58 9 

0 007668 

47 52 

0 258498 

1881 

1 13499 



34 78 

0 583882 

18 92 

1 110329 

IIQQB 

0 558693 


0 558625 


0 609221 

51 91 

0 157962 





























































51 06 

245 

0 15 

51 83 

0 158138 

41 49 

245 

0 85 

24 27 

0 909734 

22 04 

245 

HDHI 


1 159631 

34 59 

244 84 



1 561519 

36 63 

244 69 


■QQI 

0 621381 

44 


0 39 

49 43 

0 410452 


33 49 

244 69 

51 22 

244 69 

37 1 

245 


17 7 35 1 812408 




244 69 

0 62 

32 27 

0 659121 

n^gm 

245 

0 15 

52 15 

0 157791 


241 08 

1 93" 

5 44 

2 065036 

41 65 

245 

0 15 

51 9 

0 15797 


244 84 

1 46 

10 34 

1 561519 


239 98 

0 53 

37 42 

0 559935 

25 96 


0 87 

24 87 

0 921832 

39 61 


1 89 

5,42 

2 01357 


30 51 241 86 0 53 37 26 0 559042 


39 29 


242 18 


0 26 


48 23 0 271103 


58 9 


244 8^ 


0.97J125 


39 14 I 244 84 I 0.85 


27.06 

245 

0 53 

36.22 

0 558409 

41 49 

245 

0.54 

35 23 

0.571347 

41 18 

244 53 

0 16 

51 55 

0 170566 

37 88 

244 69 ^ 

0 16 ^ 

51 55 

0 170479 

32 08 1 

244 69 

2 81 

10 84 

3 005249 
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Table 7.3.5 . NSGA population for process model by 3^ factorial design on 
Moldflow Simulator for LDPE Diskcady 


— Pareto Optimal Solutions after 50 generations 1 

I Process Variables 

Responses 

Mold Tem|). 

°c 

Melt Temp. 

°c 

Inj. Time 

See 

Temp. Diff. 

Iiij. Pre. 

MPa 

59 84 

244 37 

1 24 

6 34 

13 89499 

60 

245 

0 57 

3 5 

33 73152 

58 59 

245 

1 45 

7 16 

9 658441 

58 59 

244 69 

0 85 

4 82 

24 10034 

53 57 

245 

0 28 

2 12 

45 84569 

59 84 

243 27 

0 14 

1 27 

53 21534 

58.59 

245 

0 22 

1 73 

48 55837 

60 

243 27 

1 12 

5 87 

16 95129 

59 84 

245 

0 56 

3 45 

34 17725 

53 57 

245 

0 23 

1 86 

47 97754 

49 8 

244 84 

0 85 

5 05 

24 12668 

39 76 

239 98 

0 75 

4 72 

29 45504 

53 57 

245 

0 56 

3 58 

34 16594 

59 84 

245 

1 73 

8 01 

5 864067 

57 33 

239 98 

0 22 

1.73 

50 66452 

59 84 

244 69 

1 89 

8 43 

4.66:605 

59.84 

245 

0 1 

1 08 

54 18497 

58.59 

245 

0.51 

3.24 

35 98876 

59.53 

243 9 

0.22 

1 72 

49 03319 

53 73 

243 27 

1.44 

7 26 

10 357 

59 84 

244.69 

1 89 

8 43 

4 666605 

57.02 

245 

1 72 

8 1 

6 039097 

58.59 

243 27 

0.85 

48 

24 56144 

60 

243 27 

1 45 

7 06 

10.04894 

59.84 

245 

2 47 

9 75 

5 010248 

58.27 

245 

0 12 

1 16 

53 58505 

58 9 

244 84 

1 8 

8 25 

5 264582 

60 

244 37 

1 24 

6 34 

13 89387 

60 

243 43 

1 07 

5 64 

18 41049 

58 59 

245 

1 5 

7 32 

8 89762 

59 84 

245 

0.15 

1 33 

51 90136 

58 59 

240 76 

0 85 

4 76 

25 38284 

59.84 

245 

0 12 

1 14 

53 60911 

58 27 

245 

0.37 

2 54 

41 76141 
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52 31 


59 84 


59 84 


59 53 
59 84 


60 


57 02 


57 49 


58 59 


58.59 


58 59 


59 84 


59,84 


243 9 


239 67 


244 37 
243 9 


60 

245 

59 84 

244 84 

60 

245 


224 92 


244 69 


245 


240 76 


2 


49 12126 


52 4673 


6 957197 


53 53292 


4 809429 
45 94376 


58 20898 


49 03756 


53 21534 


59 5043 1 


53 87811 


33 24719 


4 686086 


6 137174 


52 46986 


63 67374 


37.49955 


57 67107 


25 38284 


45 87907 


10 06718 


5 303533 


59 84 


59 84 


49 8 


59,84 


245 


233 24 0 85 


56 


59 06 

245 

57 02 

239 98 


243 43 


29 08792 


49 09664 


48 57517 


3 885122 


27 96214 


34 17725 


8 880527 


703 


25,01153 


7 390571 


53 57 

243 75 

0.1 

1 14 

54 6655 

53 57 

244 84 

1 62 

791 

7 368958 

59 84 

245 

0 35 

2 39 

42 78643 
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Table 7.3.6 : NSGA population for process model by 3^ factorial design on 
Moldflow Simulator for LDPE Diskcady 


Pareto Optimal Solutions 1 

I Process Variables 

Responses 

Mold Temp, 

°c 

Melt Temp. 

°-C 

Inj. Time 

_Scc 

Temp. Diff 

Ac. Inj. Time 

Sec 

60 

245 

1 51 

73 

1 61308 

60 

244 37 

001 

0 55 

0 008064 

60 

244 84 

0 01 

0 55 

0 007986 

60 

239,67 

0 01 

0 54 

0 008743 

60 

239 2 

0 06 

0,8 

0 058917 

59 37 

245 

0 01 

0.55 

0 007951 

60 

244 84 

001 

0 55 

0 007986 

59,69 

245 

0.01 

0 55 

0 007955 

58 75 

245 

0.01 

0 56 

0 007942 

60 

243 75 

0 01 

0 55 

0 008165 

49 96 

245 

001 

0 64 

0 00782 

59 84 

239 98 

001 

0 55 

0 008701 

54 98 

244 84 

0 01 

0 59 

0 007916 

54 2 

234 96 

1 51 

7 24 

1 614455 

60 

245 

0 01 

0,55 

0 007959 

49.33 

243 59 

0,01 

0 64 

0 008042 

60 

244 84 

0 01 

0 55 

0 007986 

60 

245 

02 

1 59 

0 208416 

53,1 

245 

1 51 

7 58 

1 612696 

60 

239 98 

0 39 

2 53 

0 409721 

56 86 

242 49 

0 01 

0 57 

0 008314 

60 

239 2 

001 

0 54 

0 0088 

60 

245 

0,01 

0,55 

0 007959 

59 84 

245 

0.01 

0 55 

0 007957 

60 

239 2 

0 76 

4 29 

0 811063 

49 65 

245 

0 01 

0 64 

0 007816 

60 

245 

001 

0 55 

0 007959 

59 37 

245 

001 

0 55 

0 007951 

59 84 

245 

001 

0 55 

0 007957 

59 84 

245 

001 

0 55 

0,007957 

59,84 

245 

0 01 

0 55 

0 007957 

60 

244 69 

001 

0 55 

0 008012 

60 

245 

001 

0 55 

0 007959 

49 65 

245 

0 1 

1 19 

0 10801 
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60 

245 


59 37 

239.35 


60 

244 22 



54 35 

245 

60 

239 98 


39.29 


60 


0 00801 


0 007959 
0 008773 
0 00809 


0 008703 


44 63 

239 98 

02 

1 76 

0 208892 

, 59 37 

245 

001 

0 55 

0 007951 

60 

245 

0.01 

0.55 

0.007959 


0,007672 


0.007959 


58 75 


0 003381 


0.058174 


0 007942 


59 84 


244 69 


244 69 


5 


0 007951 


0 007872 


0 00801 


0.007959 


0 


0 033119 


0 009665 


59,84 


59 84 


59.69 


58 75 


60 


39,92 


57.18 


59 69 


6 


245 


234 96 


245 


245 
244 84 


243 75 


0.007957 
0 007957 
0 007995 


0 007942 


0 033011 


0 007681 


0 007947 


0 208404 


0 007955 


0 007959 


0 007957 


0 008165 
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Table 7.3.7 : NSGA population for process model by CCD rotatable design on 
Moldflow Simulator for Polystyrene Diskcady 


Pareto-Ontimal Solutions af 

ter 50 Generations 

Process Variables 

Responses 

Mold Temp. 

°C. 

Melt Temp, 

°c 

Inj. Time 

. _ Sec 

Inj. Pre. 

MPa 

Temp. Diff, 

70 

224 06 

0 37 

25 57 

0 396703 

69 61 

225 

1 52 

10 74 

1 627452 

70 

224 69 

0 29 

26 84 

0 308767 

70 

225 

1 17 

13 25 

1 251556 

70 

224 22 

0 37 

25 47 

0 396758 

70 

224 37 

03 

26 79 

0.321252 

69 22 

224 37 

0,37 

25 63 

0.396945 

69 8 

224 37 

0.25 

27 83 

0 270901 

69 8 

225 

1 63 

10 15 

1 740106 

70 

224 69 

1 17 

13 44 

1 251446 

57,06 

225 

1 59 

14 84 

1 700359 

69 8 

219 98 

1.63 

13 08 

1 73834 

70 

224.22 

0.37 

25 47 

0.396758 

70 

224,53 

1 51 

10.92 

1 614788 

70 

224,84 

0,12 

30 31 

0.132399 

69 8 

224 37 

1.51 

1111 

1 614724 

69 61 

224 37 

1 45 

11 54 

1 552109 

69.61 

214 8 

1.91 

15 15 

2 036594 

70 

224 06 

0 37 

25 57 

0 396703 

41 57 

219 04 

0.77 

26 74 

0 820921 

70 

224 69 

03 

26 6 

0,321363 

69 8 

223 9 

02 

29 38 

0 207735 

70 

225 

1.7 

981 

1 815186 

44.9 

224 84 

0 12 

33 91 

0 132607 

69.8 

223,9 

0 95 

16 48 

1 01276 

69 8 

224.37 

02 

29 09 

0 2079 

70 

224.37 

02 

29.03 

0 20786 

69 8 

224.69 

0.37 

25 25 

0 396957 

69 8 

224 37 

1 45 

11 45 

1 552115 

70 

222 49 

081 

19 13 

0 861527 

69 8 

225 

1 63 

10 15 

1 740106 

70 

223 27 

001 

34 

0 00574 

70 

224 37 

02 

29 03 

0 20786 

69.8 

224.22 

0,14 

30 49 

0,144827 
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69 8 


63 53 


223 9 


225 


29 38 


0 018844 


1 552122 


1 012489 


0 207735 


57 45 

224 53 

70 

225 

70 

224 37 

69 8 

225 

69 8 

224,37 


70 


70 


68 24 


69 8 


70 


69 8 
69 61" 


224,37 


223 27 


223 9 


224.37 


224 22 


224 22 


6 


69 8 

224.84 

70 

224 37 

69.8 

225 

69 8 

224 37 


69 61 

224 22 

70 

224 69 


6 


224.53 


224.37 


224 37 


69 41 

225 

69 61 

225 


26 41 


26 79 


0 221774 


0 321474 


0 321252 


1 740106 


1 552115 


0 724053 


1 732715 


1 614732 


25 38 

0 396814 

34 

0 00574 


0 811929 


11 53 

1 53959 r 

33 7 

0 006541 

27 87 

0 270807 


25 57 


13 82 


0 811895 


0.396703 


1 251224 


30 66 

0 132276 

25 19 

0 396924 

10 15 

1 740106 


33 54 


32 97 


1,865159 


0 396826 
1,226361 


0 006017 


0 018797 


10 48 


3 64 


1 702537 
"l 226478 


155 





























































69 41 


69 8 


63 73 


70 


69 8 


225 


224 37 


8 


25.19 


1 22648 


0 270901 


1 614152 


0 396924 


552 


70 

204 61 

03 

36 05 

0.314142 

63 33 

225 

1 63 

1274 

1 739276 

69 8 

224.37 

1 51 

11 11 

1 614724 

70 

224 22 

0 57 

21 8 

0 610706 ' 


70 

224 69 

66 47 

225 

69 8 

224 37 

... 

69 8 


68.24 

111 18 

69,61 

225 

70 


70 

225 

69.8 

224.37 

69 61 

225 

70 

224.06 


17 87 


10 54 


7 


25 19 


0 849867 


1,764855 


0 

"0 321252 


0 396814 


0 00-M58 


0 396924 


1 552104 


28.18 0 258247 


10 81 


32 87 


26 41 


11 45 


10.39 


25,57 


1.614933 


0 018853 


0 321474 


1 552115 


1 702551 
0.396703 



156 























































Table 7.3.8 : NSGA population for process model by CCD rotatable design 
Mold flow Simulator for Polystyrene Diskcady 


Pareto Ontimal Solutions of 

ter 50 Generations 

Pi 

rocess Variables 

Responses 

Mold Temp. 

Melt Temp. 

^ 

Inj. Time 

Sec 

Temp. Diff, 

°C 

Inj, Pre. 

MPii 

70 

224 84 

0 91 

4 25 

IT i.T g 

16.28481 

70 

224 84 

1.8 

67 

9 629505 

70 

225 

0.01 

0 92 

32 96558 

70 

214 65 

0 02 

0.88 

38 29068 

70 

224 22 

0 86 

4 08 

17 28967 

70 

224 84 

1 29 

5 39 

12 28337 

70 

224 84 

0.3 

2 09 

26 5074 

70 

224 84 

03 

2 09 

26 5074 

70 

224 22 

0 91 

4 24 

16 6669 

70 

224.84 

03 

2 09 

26 5074 

57.45 

225 

0 03 

1 35 

35 11712 

70 

219.82 

1 1 

4 68 

16 97155 

70 

224.22 

0.86 

4 08 

17.28967 

70 

224.84 

1.85 

6.8 

9.551772 

70 

224.22 

0.01 

0.91 

33.44254 

70 

224.84 

0.02 

0.97 

32 77904 

70 

224.84 

0.3 

2 09 

26 5074 

70 

224.84 

0.29 

2 05 

26 74913 

70 

224 84 

0.01 

0 92 

33 06165 

69.8 

224.22 

1 43 

5.76 

1 1 69858 

70 

224 84 

0 91 

4 25 

16 28481 

69.8 

225 

0.03 

1,02 

32 46019 

69.22 

224.84 

03 

2 12 

26 7564 

70 

224 84 

0 01 

0 92 

33 06165 

70 

224 84 

0 86 

4.1 

16 90764 

70 

224.84 

0.11 

1.35 

30 57948 

70 

224 84 

0 11 

1.35 

30.57948 

70 

224.84 

0 1 

1 3 

30.84846 

70 

224.84 

03 

2 09 

26 5074 

70 

224.84 

03 

2 09 

26 5074 

70 

224,84 

0 35 

2 27 

25 55752 

70 

224.84 

0 07 

1 16 

31 66563 

70 

224,84 

0 11 

1,35 

30 57948 

70 

224.84 

0.11 

1 35 

30 57948 


on 


157 





09 


33 44254 


69 8 


57 45 


70 


70 


70 


70 


70 


70 


43,33 


66.86 


70 


70 


69.8 


70 


69,61 


70 


224 84 


224 22 


224 22 


22 


224 84 


224 84 


224 84 


219 82 


214 8 


219.82 


224 84 


224 84 


14 37241 


33 89055 


27 35293 


15 03866 


32 Oil 18 


1 29 

5 39 

12 28337 

1 29 

5 39 

12 28337 


30 57948 


12 28337 


16 28481 


31 66563 


25.55752 


26.5074 


26 69873 


30 57948 


0 11 

1,35 

30 57948 

0 11 

1 35 

30 57948 


29 90167 


0.97 

32 77904 

2.1 

26 47422 


35.67305 


30 69823 


224 84 


11 02142 
17 36281 


To 

214.8 

10 

” 224.84 


224.84 


31.95119 


25 55752 


26 5074 


31 94025 
14 96214 


158 








































































70 

225 

03 

2 1 

26 41123 

70 

224 84 

0.3 

2.09 


63 " 7"3 

224 84 

0 3 

2.29 

28 27582 

70 

224 84 


1 54 

29 52058 

70 

224 84 


1 35 

30 57948 

70 

204.76 

0 35 

1 99 

35 04617 

63 73 

224 22 

0 01 

1.08 

35 04296 

70 

223.59 

0.35 

2 25 

26.31474 

70 

224 84 



29.52058 




26.74913 
33 06165 
30 57948 



159 
































































Table 7.3.9: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 


Qptitnal Solutions after 50 Generations 


Process Variabl es 

Melt Temp. | Mold Temp. I M. Inj. Pre. I M, Inj. RatT 


Response s 

Fill Time I Cooling Time 


211.92 


206 9 


217.57 


206,9 


206 9 


208 16 


206.9 


208.16 


211 92 


208.16 


217.57 


216 94 


211 92 


206.9 


208.16 


220.71 


213.18 


201 88 


216.94 


200.94 


206.9 


201 88 


201.88 


211.92 


205.65 


220.71 


201.88 


298 04 


298.04 


294 12 


298.04 


294 12 


2 


298,04 


10000 


10000 


0.000932 


0.000916 


0 000923 


0 000927 


2,162 


2 44 


2 289 



35.18 


35 


35 


40 65 


35 


35 


35 


298 04 


298 04 



2 


172.55 


2 


298 04 


282,35 


298 04 


298.04 


298.04 


10000 


10000 


10000 


10000 


10000 


10000 


10000 


HI 


10000 


10000 


10000 


10000 


0 000914 


0.000927 


0.000927 


0 000926 


0 000927 


0.000925 


0 000926 


2 437 


2.226 


2 226 


2 243 



0 000926 

2 242 

0 001265 

2 359 


0 000919 


2 349 


0.000922 

2 289 

0 000927 

2 226 



0 000933 


0 000917 


0,000921 


0 000932 


0 000919 


0.000933 


0.000927 


0 001043 


0 000932 


0.000922 


0 000928 


0 000917 


0.001028 


0.000914 



2.304 


2 






























































































201 88 


221 96 


224 47 


201 88 


208 , 


200 63 


201 88 


201.88 


208 16 


224.78 


206 9 


201 88 


206.9 


225 73 


208 16 


217.57 


204.71 


206.9 


211.92 


208.16 


208.16 


205.65 


206,9 


201 88 


220.71 


211.92 


9701.96 




0.000941 





210.67 

35 

213 18 

35 

206.9 

35 


290 2 


298 04 


282.35 


298.04 


294.12 


10000 


10000 


9701 96 


0 000935 


0 000932 


0 000936 



35 


35 


35 


35 


35.18 


35 


35 


35 


35 


57.59 


35 


35 


35 


298.04 


298,04 


298.04 


298.04 


298 04 


298 04 


297.06 


298 04 


10000 


,96 


10000 


0.000929 


0.000922 


0,000927 


0,000927 


217.57 

35 

294.12 

209 1 

35 

1 294.12 

— 

206.9 

35 

298.04 

i 211.92 

35 

294.12 ~ 


10000 


10000 0 000927 _ 


10000 0.000925 


298.04 


297.06 


298 04 


0 0.000932 


10000 1 0 000917 


2 226 


2 162 


2 243 


2 146 


2 162 


2 21 


2 162 


2 242 


2 44 


2 242 


2,357 


2198 


2.226 


2.289 


2 227 


2 226 






2 226 


2 289 


2.242 


221 


2 993 


2 162 


2 393 


2 289 


161 





































































































208 . 


206 9 


206 9 


217 57 


225 73 


208 16 



298 04 


298 04 


10000 


10000 


10000 


0 000915 


0 000925 


0 000929 


2 242 


2 194 



162 


























































Table 7.3.10: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 


Pareto Optimal Solutions after 50 

Generations 

Process Variables 

Responses 



279 69 

76 47 

279 69 

79 82 

280 

79.12 


280 


280 


279 69 


280 


280 


280 


280 


28 


280 


280 


280 


280 


28 


280 


75 94 


74.18 


76 47 


79 12 


79 29 


73 47 


79.12 


79.29 


a 


77 88 


79 12 


8 


79 29 


79 29 




77 71 


79 29 


55 88 


144 12 


114 71 


9888 24 


9552.94 


9441.18 


206 86 


105 88 


178 43 


245 1 


144 12 


71 57 


148.04 


9925.49 


9888 24 


9441.18 



2 


245 1 


58 82 


175 49 


50 


9 

Y' 


86 27 




9552.94 


9925 49 


9925 49 


9441 18 


9925.49 


9552 94 


9776 47 


9515 69 


70.273 


71 777 


62 974 


62 498 


57 83 


75 938 


65 241 


72 987 


pBE 

SI 


70 646 


8360 78 


9776 47 


9552 94 


9739 22 


54 727 


73 387 


0 00246 


9888 24 

75 881 

0 000983 

9813 73 

57 159 

0 001852 

9552 94 

73 484 

0 001072 

9888 24 

70 143 

0 001211 

9552 94 

66.19 

0 001414 


0 001212 


0 00115 


0 002595 


0 001603 


0.001816 


0.000981 


0 00144 


0 001097 


0 . 

000933 

0 

001238 

0 

001702 

0 

001189 

0 

001603 

0 

001165 



0 000983 


0 00095 



0 002182 


0 001081 


035 





























































































279 69 


280 


280 


28 


280 


279 69 


280 


280 


280 


280 


28 


280 


280 


28 


280 


279 69 


280 


280 


280 


28 


279 69 


279 69 


277.49 


277 18 


280 


77 88 


79 12 


146 08 



79 82 


79 12 


80 


79 82 


8 63 


50 98 


208 82 


9925 49 


9925 49 


9441 18 


9105 88 


10000 


9552 94 


7466 67 


76 305 


77 105 


0 00123 


0 002124 


0 001418 


0 001702 


0 001079 


0 000956 


73 74 

0 001087 j 


0 00 

lEa 

66 628 

0 00 

1387 

51669 

0 00246 


0 000983 




.2 


55 88 


9552 94 

69 731 

9515 69 

62 559 

9925 49 

57 36 


0 001593 


0.001845 


280 

79 82 

144.12 


0 001438 

280 

79 82 

58.82 

9925.49 

57 721 

0001813 

279 69 

79 82 

58 82 

7466 67 

54 97 

0 00238 

277 49 

79 12 

245 1 

9627 45 

78 22 

0 000943 

280 

- - I 

80 

147 06 

9515 69 

70 004 

0 001227 


79 12 


80 


245 1 


7 06 


9552 94 


77 32 


70.004 



76.47 


76 29 


79.29 


79 29 


77 88 


2 


94 12 


208 82 


208 82 


90 2 


50 


55.88 


224 51 


182.35 


9739 22 


10000 


10000 


9925.49 


8360 78 


8994.12 


63 379 


0 001227 


0 001383 


0 001538 


0 000975 




9515 69 



7 

162.75 

9552 94 




7 


79 29 


73 47 


62 75 


206 86 


9739.22 


9478.43 


58 716 


0.001551 


0.002182 


0 001977 


0 000954 


0 001076 


0 001151 


0 00181 



79 12 

206 86 

9552.94 


79 12 


9552 94 





wmmi 

m 

0 001002 

1 

0.001898 


9552.94 


9552 94 


66.19 


69 731 


0 001817 


0.001007 


0 001414 


0 001238 


164 











































































































































280 

73.65 

62 75 

10000 

58 567 

0 001774 

280 

76 29 

208 82 

9739 22 

75 856 

0 000986 

280 

80 

90 2 

9515 69 

62 559 

0 001593 

280 

79 82 

118 63 

8696 08 


280 

79 12 

208 82 


0 000983 


0 001047 


269 96 

77 88 

146 08 

9925 49 

72 816 

0 001227 

280 

79 12 

245 1 


77 321 

0 000948 

280 

80 

146 08 


0 002971 

280 

79 12 

92 16 

9515 69 

62 873 

0 001579 

279 69 


177 45 1 

10000 

73 916 

0 001063 



9105 88 

72 351 

0 00114 




8994 12 

72 38 

0 001158 

278 75 

80 

177.45 

10000 

74 002 

0 001063 


279 69 


280 


279.69 


239.53 


279 69 


280 


74 18 


2 


79 12 


79 29 


91 18 


55 88 


4784 31 


9552 94 


9925 49 


9925 49 


79 82 


63 001 


81 585 


57 481 


1 280 

79 82 

83 33 

280 

79 12 

208 82 

279 69 

79 82 

58 82 

280 

79 12 

191 18 

280 

79 12 

177 45 


9515.69 


9515 69 


9850 98 


9776 47 


9441 18 


6 


75 955 


57.675 


57 


72.912 


0 003427 


0 001883 


0 001553 
0 001305 


0.001862 


0 001425 


0 001648 


0 00183 


0.001047 


0001101 


165 




















































































Table 7.3.11: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 


Pareto Optimal Solutions after 50 Generations 


Responses 


256 


256 . 


256 16 


257 1 


257 


0 000908 


259 92 


259 92 


257 1 


259 61 


257 1 


256 16 


257 


259 61 


257 , 


257 1 


256.16 


259 92 


256 16 


37 82 


8 


35 


257 1 

35 

259 92 

35 

257 1 

35 


40 65 


35 


35 


35 


35 


35 


270 59 


274 51 


270 59 


270 59 


.59 


270 59 


270 59 


270 59 


270 59 


270 59 


270 59 


256 16 

35 

257.1 

35 


207 84 


59 


259 92 

35 

27 

259 9 

'2 

35 

27 


6 

35 

27 


27 


278 43 


166 


10000 


10000 


10000 


10000 


10000 


0000 

0000 


10000 


10000 


10000 


10000 


10000 


10000 


0 115 

0 000908 

0 127 

0 00091 

0 114 

0 000908 

0.114 

0 000908 

0.114 

0 000908 

0 111 

0 000909 

0 114 

0 000908 

0 111 

0 000909 


270 59 

10000 

0 114 

0 000908 

270 59 

10000 

0 111 

0 000909 

270 59 

10000 

0 132 

0 000912 

270 59 

10000 

0 111 

0 000909 

145 1 

10000 

0 18 

0 001268 


0 117 


0 114 


0 115 


0 114 


0 115 


0 115 


0 12 


0 115 


0 000909 


0 000908 


0 000908 


0 000915 


0 000908 


0 000909 



0 000908 


0 000908 


0.000908 


0 000908 


0.111 

0 000909 

0 115 

0 000908 



0 000909 


0 000909 

5 

0 000908 


0 000909 

) 

0 00091 







































































































257 41 


257 1 


257 


257 1 


259 92 


259,92 


257 41 


257.41 


259 92 


259 61 


259 92 


257 


256 


257 1 


259 92 


258 98 


259.92 


259.6 


259.92 


259 92 


259 92 


256 16 


35 18 


35 


35 


35 


35 


35 


35 


270 59 
2”76 59 


10000 

■im 

0 000908 

8807 84 

■B 

0 000999 

10000 

0 113 

0 000908 

10000 

0111 

0 000909 

10000 

0 133 

0 000908 


0 113 

0 000908 


35 

27 

35 

27 


270 59 


239 22 


59 


27 


59 


.59 


270 59 


270 59 


10000 


10000 


10000 


10000 


10000 


10000 


10000 


59 


270 59 


270.59 


270 59 
270 59 
270 59 


10000 


0 


5231 37 


10000 


10000 


10000 

10000 


0 113 


0,114 


0 114 


0 111 


0 113 


5 


0 111 



8 

27 

EBB 




0 111 


0 114 


6 


15 


0 373 


0 111 


0 111 


0 111 


0 000908 


0 113 

0 000908 

0.11 

0 000927 


0 000908 


0 00091 


0 000909 


0,000908 


0,000908 


0.000916 


0 000909 


.59 

10000 

0 111 

0.000909 


0 000909 


0.000909 


0 000908 


0 000908 


0 000908 


0 000909 


0 000908 


0 002217 


0 000909 


0 000909 


0 000909 


0 000909 



35 

270 59 

10000 

0 115 

0 000908 


35 

270 59 

10000 

0.113 

0 000908 

259 92 

35 

274 51 

10000 

0 111 

0 00091 


167 












































































































259 61 

35 18 


259 92 

35 


255 84 

35 

270 59 

257 41 

35 

1^^ 

259 92 

35 

270 59 


35 

270 59 

249 88 

35 

270 59 

259 92 

35 

270 59 

257 1 

35 

270 59 

259 92 

35 

270 59 

258 98 

35 " 

Kmyili 

257 41 

35 

270 59 

257 1 

35 


256.16 

35 


259 92 

40 65 

270.59 

259 92 



259 92 

35 

^^1 

256 16 

35 

270 59 

259 92 

35 

270 59 

257 41 

35 

270.59 

256 16 

35 

270.59 

^agnail 


WSmml 

257 1 

35 


259 92 

35 

270 59 

257 1 

35 

270.59 

256 16 

35 

270.59 


10000 


0 112 


0 000909 


10000 


0 111 


0 000909 


10000 


0 116 


0 000908 


8807.84 


0 141 


0 000999 


10000 


0.111 


0 000909 


10000 


0.114 


0 000908 


10000 


0 127 


0 000908 


10000 


0 111 


0.000909 


10000 


Bililrfijgli?;! 

10000 

KUIB 

0.000909 



0.000909 

10000 

■UH 

0 000908 

8807 84 

Igf 

0 000999 

10000 

0 115 

0 000908 

10000 

0 132 

0 000912 

10000 

0 111 

0 000909 

10000 

0.111 

0.000909 

10000 

0 115 

0 000908 

10000 

0.111 

0 000909 

10000 

0 113 

0 000908 


0.12 

0000916 

10000 

0 1 15 

0 000908 

10000 

0.114 

0 000908 


mam 

0 000909 

mil 


0 000908 

■EElBil 

0 115 

0 000908 


168 































































































Table 7.3.12: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 


Pareto Optimal Solutions after 50 Generations 


Process Variables Responses 



200 


200 627 


2 


200 


2 


215.059 


215 059 


245 176 


215 686 


240 157 


200 


200 


200 


245 804 




35 529 



51 961 

574 51 

52 325 

2 191 

50 

n 


51 6 

2 2 


574 5 
574.5 


574.51 


574 51 


51 956 




48.068 


255 216 

35 529 

215 059 

35 529 

201 882 

35 529 

200 

35 353 

255 843 

35 529 


51 961 


51 961 


51 961 






574 51 


574 51 


574 51 




51 752 


52 327 


37 691 




35 529 


35.529 


35,529 


51.961 


67.647 


51 961 


500 


574.51 


574.51 


574 51 


574 51 


43,697 


52 325 


51 956 


51 956 


47 383 



872 549 


49 896 


240 157 

35 529 

215 059 

35 529 

240 784 

35 529 


2 191 






574 51 

47 548 

2 377 



2215 


2 191 


2 799 


2 257 


77 


2,36 


2 377 


2,7 


2 385 


2 648 


2 191 


2.191 


2 191 



2 377 


loom 

62 355 ' 

2 182 


39 49 

2701 

723 529 

48 204 

2 376 


2 374 



574 51 


51 755 


2 203 

























































































200 

35 529 

50 

723 529 


200 


275 922 


205 02 


255 843 


252 078 


210 039 


205 02 


36 235 


35 529 


35 529 


35 529 


35 529 


35 529 


51.961 


67 647 


51.961 


5 


51 961 


52.609 

2 191 

42.897 

2.788 

51 946 

2 195 

33 367 

2 965 


537 255 


52 159 


2 255 


2 8 


2 822 


2317 


2 255 


2 191 


200 

35 529 

50 

200 

35.529 

50 


35 529 

51 961 


251 451 


215 059 


205 02 


200 627 


41 176 


35 529 


35 529 


96 


50 


50 


50 

181 373 
50 


200 

35 529 

50 

240 784 

35.529 

50 

205 02 

35 529 

50 


55 882 


574 51 


574 51 


723 529 


574 51 


872 549 
872 549 
537 255 


1133 333 


57.065 


37 828 


52 327 


39011 


47 548 


48 152 
69 941 
51.6 


50 447 


38 4 


40 852 


2 


2 184 


2 794 


2 191 


2 828 


2 377 


2.405 
2 236 
2 2 


2 


2.792 


2.669 


2 65 


2651 


200 


200 3 


35 529 


50 


52 941 


574 51 


574 51 


574 51 


37 538 


51 956 


52 412 


2 794 


2 191 


2 195 


200 627 


2 


215 059 


200 627 


200 


35 529 


35 353 


50 


51 961 


537 255 


5343 137 68 657 


2 196 


2 165 
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Table 7.3.13: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 


Pareto Optimal Solutions 50 Generations 


Process Variable Responses 



261.49 

35 

143 14 

256 47 

35 

268.63 

241 73 

35 

268.63 



243.92 


241.41 


244.55 


48 


261 49 


248 


268.63 


252.94 
268 63 


268 63 


268.63 


0 


10000 


0 


QB 

EH 



264.71 


268.63 


205 88 


10000 


10000 


10000 



0 142 


62 


0 115 


0 149 


0 115 


0.12 


0 142 


5 


0 14 


5 


0 114 


0 268 


0 179 


0 133 


0 277 


268 63 


268 63 


268 63 


10000 


10000 


8807.84 



3 

35 

259 8 

245 4 

9 

35 

260 78 


0 131 


0 149 


0 152 


0 14 


2814 


2 765 



2 742 


2 645 



2,342 


2 646 


2 685 


2 346 


2 648 


2 685 

2 sT 

2 685 
2 623 





























































246 43 

35 

268,63 

10000 

0 135 

2 67 

253 96 

35 

268 63 

10000 

0 119 

2 742 

mm 

35 


10000 

0 135 

2 67 


256 47 


246 75 

35 

243 92 

35 

241 73 

35 

241 73 

35 


253 96 
251 45 
216 63 


266 82 
253 96 
256 47 
246 43 


35 

35 

40 65 


35 

35 

35 

35 


268,63 


268,63 


263 73 


268 63 


268 63 


268 63 
268 63 
268 63 


268 63 
268 63 
268 63 
264 71 


ggggi 

mm 

mm 


10000 


10000 


10000 


10000 

10000 

10000 


0 


0 134 


0 144 


0 149 


9 


0 119 
0 123 


0 288 


0 106 


0 119 
0 115 
0 136 


2,673 


2 645 


2 623 
2 623 


2 742 


2718 


2,408 


2 856 


2 742 
2 765 " 

’YoT' 


256 47 

35 

267 65 

10000 

0 114 

2 765 


35 

268 63 

9701 96 

0 141 

2,67 

251 45 

35 

268 63 

10000 

0 123 

2 718 


35 

||g^Qg|[| 

10000 

0 115 

2 765 


35 


10000 

0 107 " 

2 837 


35 

268 63 

10000 

0 135 " 

2 67 


268 63 


10000 



174 
















































Table 7.3,14: NSGA population for process model by Box - Behenkin design 
on CMOLD simulator for polypropylene Tape casing 
) Pareto Optimal Solutions 50 Generations I 


Process Varia bles 

Melt Temp, Mold Temp. M. Inj. Pre. M. Inj. Rate 


Responses 
Inj. Pre. I Temp Diff. 


280 
280 
280 
280 
280 
274 98 


37 12 
37 12 

36 06 

37 12 
35 71 
37 12 

-7 


MPa 
53 92 

' 50 
50 
50 
50 
50 
50 


53 92 


50 


Cm3/scc 
6162 75 

5566 67 
1990 2 
5566 67 
7354 9 
5566 67 
5566 67 


MPa 
51 99 

49 98 
38 61 

49 98 
53,65 

50 86 
49 98 




35 88 

50 

37 65 

57 84 



2623 53 


5566.67 


5305.88 


42 11 


49,97 


49 33 


6162 75 

66 15 

0 23046 

5566 67 

49 98 

0381816 

5566 67 

49 98 

0 383219 

7354 9 

53 64 

0 269891 

5566 67 

52 99 

0 352529 

6162 75 

51 99 

0 332573 



0 332573 
0 383219 
0 784328 
0 383219" 
~0 262789 
0 42146 
0 384605 


0 690159 


0.403297 


0 406234 


0 383219 


0 326856 


0381816 


0 200276 


0 28916 


0 187216 



0269891 


0 378958 


55 97 

0 213825 

51 35 

0 332749 



0 262613 


0 384605 


7950 98 

55 92 

0 23306 

5566 67 

51 22 

0 37298 
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37 12 

50 

6013 73 

51 02 

0 349043 

IIBB 

37 65 

50 



0 385292 

280 


50 

9143 14 



280 

BOB 

53 92 

6162 75 

51 99 

0 326856 



50 

5566 67 

58 14 

0 770847 

280 

37 47 

50 

5566 67 

49 98 

0 384605 

267 45 

37 47 

50 

7354 9 

55 89 

0 358446 

280 

37 12 

50 

5566 67 

49 98 

0 383219 

280 

36 06 

50 

5566 67 

49 98 

0 378958 

280 

35 71 



50 6 

037161 

280 




IQQII 

0 695206 



50 

1692 16 

HBEBI 

0 832959 

HQIQHl 

37 29 

51 96 

6460.78 

52 3 

0 316261 

280 

37.29 

53 92 

7354 9 

54 31 

0 263308 

280 


50 

5268 63 " 

49.23 

0 424107 

280 

35.88 

50 

5566.67 

49.98 

0 378233 

IIIIIIIQI^II^ 


■lll^lll^ 






■SraMI 



0 258352 

280 

37.12 

50 

5566 67 

49 98 

0 381219 

280 

37 12 

50 

5566 67 

49 98 

0 383219 

heibih 


50 

■IBB 


0 834345 


36 06 

50 

6758 82 

52 57 

0 295966 






0 40296 

269 96 

■aai 

II^QIIIIIIII 

9143 14 

Sill 

0 27171 

280 

WEIM3M 


2623 53 

42 11 

0 690159 

280 

yiM 

53 92 

8845 1 

56 38 

0 213204 
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confirmed that such solutions may be verified subsequently by comparing the “optimum” 
given by the model with a subsequent run of the simulator 

3 Product quality in injection molding is a function of the material used, the mold design, 
cooling system design and processing conditions The processing conditions that are 
important for quality of the injection molded part are mold temperature, melt temperature 
and fill time Quantitative measures of the injection molded part quality are developed 
These quantitative measures constitutes the objective function that must be minimised The 
DOE framework can isolate the sigmficant process factors from those that are not 

4 Second-order regression models appear to reasonably good to represent the process 
models for flow simulation outputs from a injection molding simulation software 
However, the models are not perfect and require verification runs before the final results 
may be accepted 

5 Population based “genetic algorithms” is well suited to injection molding process 
parameter optimisation Since the GA gives a group of solutions, decision maker may 
subsequently choose a particular solution depending on his/her requirements and 
preference 

6 Multiple objective process optimisation using GA has been done in this work for injection 
molding The concept of Pareto-optimality was used to solve the Bi-objective problems 
using “non-dominated sorting genetic algorithm” (NSGA) To the best of the 
investigator’s knowledge, this has not been attempted before (see for instance, Kumar A, 
1998) 
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Scope for future work 


• The present work has used sample CAD models for its studies CAD models of injection 
molding parts from a real industry and the optimisation frame work developed in the 
present work for process parameter optimisation should be combined and tested more 
extensively 

• In this study we have used a second regression models One may go for the higher order 
models for the better results 

• In this study we have done bi-objective optimisation of the injection molding process 
parameters One may do the multi-optimisation of the injection molding process 
parameters 

• In our study we were not able to evaluate the effect of NSGA parameters like on the 
distribution of population on different fronts This points towards the scope for a detail 
study of NSGA parameters 

• The quality of the injection molded part is affected by the mold design and cooling system 
design apart from the process conditions. Therefore, one may go for the optimisation of 
the parameters of injection molded design and cooling system design using GA. 

• An enhanced mult-objective GA (ENGA) (Srinivas, 1 998) which is proved to be faster in 
convergence than NSGA, ENGA can be used for the multi-objective optimisation 
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